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Abstract 

We present methods for evaluating the Feynman parameter integrals associated with the pen- 
tagon diagram in 4 — 2e dimensions, along with explicit results for the integrals with all masses 
vanishing or with one non-vanishing external mass. The scalar pentagon integral can be expressed 
as a linear combination of box integrals, up to 0(e) corrections, a result which is the dimensionally- 
regulated version of a D = 4 result of Melrose, and of van Neerven and Vermaseren. We obtain and 
solve differential equations for various dimensionally-regulated box integrals with massless internal 
lines, which appear in one-loop n-point calculations in QCD. We give a procedure for constructing 
the tensor pentagon integrals needed in gauge theory, again through 0(e). 
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1. Introduction 

The search for new physics at current and future hadron colliders demands that we first refine 
our understanding of events originating in known physics, most importantly in QCD and QCD- 
associated processes. To date, the matrix elements for all pure QCD processes with up to seven 
external legs are known exactly at tree-level [1] allowing the computation of events with up to five 
jets in the final state. (Various techniques [2] allow one to approximate cross sections with more 
jets.) Because the perturbation expansion for jet physics in QCD is not an expansion strictly in 
the coupling constant, but is rather an expansion in the coupling constant times various infrared 
logarithms, radiative corrections play an important role in matching theoretical expectations to 
experimental data. The calculation of radiative corrections requires of course the computation of 
loop corrections to the basic tree-level partonic processes. Thus far, the one-loop corrections are 
known only for the most basic processes, matrix elements with four external partons [3]. To go 
beyond these basic processes in the computation of radiative corrections in pure QCD (for example, 
to calculate the next-to-leading order corrections to three-jet production at hadron colliders), one 
must calculate five-point one-loop amplitudes in a theory with massless particles; and these in turn 
require the computation of one-loop Feynman parameter integrals with five external legs, within 
the dimensional regularization method. To discuss one- loop corrections to five-point amplitudes 
with external W and Z bosons, at least one of the external legs must be massive. In the present 
paper we address the computation of such dimensionally-regulated pentagon (and higher-point) 
integrals. Recently the techniques described in this paper have been used in the calculation of the 
one-loop helicity amplitudes for five external gluons [4]. 

Various authors [5,6,7,8] have discussed the computation of pentagon integrals that can be 
carried out in dimension D = 4 (i.e. that have neither soft nor collinear infrared divergences). 
In particular, Melrose [5] and independently van Neerven and Vermaseren [7] were able to express 
pentagon integrals as linear combinations of five different loop integrals with four external legs. Such 
box integrals (which, with external masses but no internal masses, are also required in radiative 
calculations in QCD) can be calculated readily in dimensional regularization, by direct integration 
or in terms of hypergeometric functions, if the number of masses is not too large. 

The techniques of Melrose and of van Neerven and Vermaseren do not apply directly to 
dimensionally-regulated integrals, however, and the required pentagon integrals have not yet been 
presented in a closed and useful form, which is to say with all poles in e = (4 — D)/2 manifest, and 
with all functions of the kinematic invariants expressed in terms of logarithms and polylogarithms.* 
Here we will provide such an expression for the basic scalar integral. We employ a set of equations 

* We have been informed that R. K. Ellis, W. T. Giele, and E. Yehudai [9] have recently evaluated the pentagon 
integrals by an independent technique. 
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derived in a separate paper [10]. These equations actually apply more generally to dimensionally- 
regulated one-loop n-point integrals; they can be used as a starting point for the reduction of an 
(n > 5)-point integral to a linear combination of boxes. For the pentagon integral (n = 5), the 
equations are the dimensionally-regulated analogs of equations derived in references [5,7]. In this 
paper we will use the equations to obtain explicit expressions for the pentagon with all lines mass- 
less, and for the pentagon with one massive external line, up to 0(e) corrections. Such integrals 
are of use in the calculation of next-to-leading-order contributions to processes such as gg — > ggg 
and Z -> qqgg. 

Besides the scalar pentagon integral, in QCD one requires tensor integrals — loop integrals 
with up to five powers of the loop momentum inserted. In the string-based technique [11,4] for 
evaluating QCD amplitudes, one obtains directly integrals over Feynman parameters rather than 
loop momenta. Tensor integrals correspond in this framework to the insertion of polynomials in the 
Feynman parameters into the numerator of the integrand. In order to construct an integral table 
that meshes well with this technique, we choose to work in terms of the Feynman-parametrized 
integrals. This approach also lets us take advantage of an observation that appropriate derivatives 
of the scalar pentagon insert Feynman parameters into the numerator of the integrand. Thus the 
scalar pentagon may be used as a generating function for all the tensor integrals. 

In the more usual momentum-space approach to tensor integrals, one performs a Brown- 
Feynman [12] or Passarino-Veltman [13] reduction, solving a system of algebraic equations for the 
tensor integrals. For example, integrals with just one loop-momentum inserted in the numerator 
are reduced to a linear combination of scalar integrals [14]. The counterparts of these equations 
exist for Feynman parameter integrals. In particular, integrals with just one Feynman parameter 
inserted in the numerator can be expressed as a linear combination of scalar integrals. If one 
now equates these expressions to the above-mentioned derivative representations of the same one- 
parameter tensor integrals, one obtains a set of first-order partial differential equations for the 
scalar integral. Thus an alternate approach to determining the scalar pentagon is to solve a set of 
differential equations. The differential equations are also an efficient way to obtain various infrared 
divergent scalar box integrals, with massless internal lines but with 0, 1, 2 or 3 external masses. 
(Most of these box integrals have been obtained previously by other techniques.) Together with 
the infrared finite box integral with four external masses [6] , for which a compact form has recently 
been provided by Denner, Nierste, and Scharf [15], these constitute the set of box integrals required 
for computing one-loop ra-point amplitudes in QCD without quark masses, for any n. (These box 
integrals will appear both in the recursive determination of higher-point diagrams [5,7,10], and as 
diagrams in their own right.) 
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The partial differential equation approach just described is reminiscent of similar procedures 
for performing two-loop and higher-loop integrals (usually with fewer external legs) [16]. However, 
the latter manipulations have generally been carried out in terms of either a momentum-space 
or a configuration-space representation of the integrals, in contrast to the Feynman parameter 
representation used here. 

The rest of the paper is organized as follows: in section 2, we introduce the Feynman- 
parametrized n-point integrals, in particular the pentagon and box integrals, and we make a change 
of integration variables and kinematic variables that allows the tensor integrals to be expressed as 
derivatives of the basic scalar integral. In section 3 we present an alternative derivation of the set of 
algebraic equations derived in ref. [10]. One of these equations can be used to determine the general 
n-point scalar one-loop integral recursively, as a linear combination of (n — l)-point integrals. (For 
n > 7 there are some subtleties, as explained in appendix VI.) The other two equations are useful in 
the calculation of tensor integrals, given the scalar integral. Also, in combination with the results 
of section 2 they give partial differential equations for the scalar integral. In section 4 we begin by 
illustrating the general derivation of the partial differential equation in section 3, using the simple 
example of the box integral with all massless external legs. We then solve the differential equations 
for box integrals with 0, 1, 2 or 3 massive external legs. In section 5 we use one of the algebraic 
equations derived in section 3 to obtain explicit formulae for the pentagon with all massless exter- 
nal legs, and with one massive external leg. In section 6, we describe how to obtain the (tensor) 
pentagon integrals with Feynman parameters in the numerator, through O(e ). 

For the reader's convenience, we have collected our results for the scalar box integrals and 
for the scalar and tensor massless pentagon integrals in appendix I. In appendix II we show that 
when the integrals are infrared finite, our results for the scalar pentagon integral reduce to the 
non dimensionally-regulated result of van Neerven and Vermaseren [7]. Appendix III presents an 
argument (verifying an observation of Ellis, Giele and Yehudai) which shows that the approach 
of section 6 generates all tensor pentagon integrals needed in gauge theory calculations. In ap- 
pendix IV, we compute an integration constant for two- and three-mass boxes. In appendix V, as 
another illustration of the partial differential equation technique, we obtain a manifestly symmetric 
expression for the triangle integral with all three external legs massive, to all orders in e. (To O(e ), 
such a formula has been obtained in ref. [17].) In appendix VI, we discuss subtleties that arise in 
obtaining scalar integrals for n > 7, and in appendix VII, we derive and discuss formulas for tensor 
integrals for both the pentagon and hexagon diagrams. 
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2. Properties of Feynman Parameter Integrals 

In this section, we shall show that Feynman parameter integrals with Feynman parameters 
inserted in the numerator of the integrand (which arise from tensor integrals) are given by appro- 
priate derivatives of the basic scalar integral. We present the particular cases of the massless box 
and pentagon integrals in more detail. 

For convenience, we assume here that the masses for all internal lines vanish. (The extension 
to nonvanishing internal masses is entirely straightforward [10].) Then the n-point scalar one-loop 
integral in 4 — 2e dimensions is 

_ 2e r i 



, (2vr) 4 - 2£ p 2 (p - h) 2 (p -k 1 - k 2 ) 2 ■■■(p-k 1 -k 2 k n -i) 2 

where hi, i = l,...,n, are the external momenta and n is the usual dimensional regularization 
scale parameter. Performing the usual Feynman parametrization, and integrating out the loop 
momentum, we obtain 

/„ [1] = r(n - 2 + e) J 1 d n ai 5(1 - £^)— , (2.2) 

' J2i,j=i Sijdidj — is 

where 

I n [l] = (-l) n+1 i (^f- e ^l n (2.3) 
is the basic n-point parameter integral, the symmetric matrix Sij is defined by 

s ij = ~^(ki-\ \-kj-x) 2 , i^j; Su = 0; (i, j are mod n); (2.4) 

and where we have put in the is explicitly. The poles in I n produced by the T function prefactor 
are ultraviolet ones; the remaining poles represent infrared divergences. In explicit calculations 
of cross-sections, they will ultimately cancel corresponding poles arising from soft and collinear 
emission of particles in (n + l)-point tree- level processes. 

We shall use the notation I n [P({oj})] to denote an integral in which the polynomial P appears 
in the numerator of the integrand, 

I n \P({a t })) = T(n - 2 + e) J* d n a % 5(1 - n - 2+£ • (2-5) 

^2i,j=i Sijdidj — ie 

In QCD calculations, one encounters integrals of this form, where the degree of P is less than or 
equal to n. 

For the box (four-point) integral, the "scalar denominator" is 

4 

Sijdidj = — sa\a^ — ta 2 ai — m\a\a 2 — m^a 2 a^ — m^a-ja^ — m^a^ai , (2-6) 



where s = (k x + fc 2 ) 2 and t = (k 2 + k 3 ) 2 , and m 2 are the masses of the external legs (some or all of 
which may vanish). For the all-massless pentagon integrals, 
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^ Sijdidj = -£120103 - 5230204 - ^340305 - s 45 a 4 ai - s 51 a 5 a 2 , (2.7) 

where s^j+i = + k i+ i) 2 . A nonzero mass for external leg 5 would add a term —m^a^a-i 
to (2.7). All external indices are understood to be taken mod n for the n-point function. We 
will present our results for kinematics in the Euclidean region, where all momentum invariants Sij, 
m 2 are negative. In this region, the scalar denominator is always positive, and the integrals are 
purely real, which simplifies the resulting expressions. We define the integral for physical values 
by analytic continuation from the Euclidean region; the analytic continuation back to the physical 
region should be understood implicitly in all formulae presented below, and we shall henceforth 
leave the ie implicit. 

Following 't Hooft and Veltman [6], we make the change of integration variables in (2.2), 

UiUi 



a n = 



3 = 1 a 3 U 3 



no sum on 1, 



■13 
a n ( 1 



(2.8) 



En 
j=l a 3 U 3 

Assuming that all on are real and positive, the integral becomes 



n-4+2e 



/■i *(i-E«i) (n; =1 «,) feu 

This form for the integral is most useful if we can also define the on in such a way that all of 
the dependence on the ctj-variables is scaled out of the denominator. Let us define the a,, and 
simultaneously a matrix p, through 

Sij = • (2.10) 

aiOij 

The elements of the matrix p io - are to be thought of as additional kinematic variables, independent 
of the aj. (In specific cases many of the elements pij may be taken to be pure numbers.) 
For the four-point integral described by the denominator (2.6), we can choose 

1 1 2 m? 2 ^2 2 ^3 2 ^4 

s = , t = , m l = , m 2 = , 777,3 = > m 4 = 

a.1013 a2tt4 a\oi2 ' OL2(x-i ' ot^a^ a^a\ 

(2.11) 

(Other choices are also possible; see section 4.) Equations (2.11) do not have a unique solution in 
terms of the cv One simple solution is a,\ = a 3 = a 2 = «4 = However, we would 
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like all four cnj variables to be independent of each other, so that we can use the scalar integral as 
a generating function for integrals with insertions of all four Feynman parameters a, . Therefore we 
consider the at to be general solutions to equations (2.11), with no other constraints on them. The 
set of independent kinematic variables corresponding to the choice (2.11) is then {a^rhf}. 
For the all-massless pentagon integral, the unique solution to 

1 



*>i,i+l — 

is 



(2.12) 



S23S34 / S34S45 / S45S51 

oci = \ , a 2 = \ , «3 



S45S51S12 V • S 51«12S23 V • S 12«23S34 , , 

(2.L6) 



S51S12 / S12S23 

"4 = \ , a 5 = J . 

S23S34S45 V S34S45S51 

Because we have taken the to be negative, the on are real. No additional kinematic variables 
are necessary for the massless pentagon. 

With these choices of c^, the four and five point scalar integrals become 

2e 



»(n-)r 



/ 4 [i] = r(2 + e ) K / d*u 



[u\Us + U2U4 + m\uiU2 + 771^2^3 + I7I3U3U4 + rh\u4U\Y 



h[l] = r(3 + e)^f[a/) jf d 5 u 



2„, n. _L™2„, „, _L™2„, „, _!_ ™2„, „, l 2 + £ 

/./II I / ; -i III / 111 — |— I / / ' ' 11.11, . 

' 1 2. 



[uiti 3 + U2U4 + n 3 ti 5 + u 4 ni + u 5 u 2 ] 3+e 

(2.14) 

Further examples of the {a^pij} change of variables are to be found in sections 4,5 and 
appendix V. 

It will be helpful to define the reduced integrals 

In[P({ai})} = (n«i) In[P({ai/ai})] . (2.15) 

As we will see, dividing out the factors of on connects the tensor integrals more simply to the 
scalar integral. For the scalar integrals (that is when the polynomial is simply 1), we will use the 
abbreviated notation I n = I n [1] . 

We can use the 't Hooft-Veltman form of the tensor integrals (2.15) to obtain derivative rela- 
tions for them. Let P m ({di}) denote a homogeneous polynomial of degree m. Then for the massless 
box integral, the change of variables (2.8) gives 

/ 4 \ - m + 2e 

f 1 . Pm{{Ui}) [Y jj = i a 3 U j) 

U [P,MoM = r(2 + e) I A, 6(1 - E,«i) lui { U3 + u 2 u i?+ ' ' (2 - 16) 



which we can rewrite in terms of derivatives acting on I4, 



h r„,t{«,}). = r( |.,, , ".;, 2 ' ; r,J{^-\ /,[!] (2.1.7) 



r(l + 2e) m Vl 5 « 
Similarly, the change of variables (2.8) leads in the five-point case to 



f 1 . P m {{ui}) U2j=i a i u j , 
h[P m ({ai})} = r(3 + e) / d 5 U 8(l-^ iUi ) ^- , (2.18) 

JO [U1U3 + U 2U4 + U 3 U 5 + U4U1 + U 5 U 2 \ 

which we can write as follows, 



l-m+2e 

3+? 



h[P m {{oi})] = r( L m o + f €) ^f(^M / B [l] • (2.19) 



T(2 + 2e) \1 9q! 

These equations hold when there are external masses as well, provided that one holds fixed the 
matrix p defined in (2.10) when differentiating with respect to c^. The result for the general n-point 
integral is 

U lPMai})] = r( " r - 3 _-'; 2 + e) 2e) P„ ({£}) /. [1] ■ (2.20) 

Equation (2.20) allows one to obtain tensor integrals by different ating the basic scalar integral. 
Certain subtleties do arise in this approach; they will be dealt with in section 6 and in appendix VII. 

Using equations such as (2.17), (2.19), and (2.20), one can translate an algebraic system of 
equations for integrals with Feynman parameters inserted, into a system of partial differential 
equations for the basic scalar integral; in principle one can then solve the equations for the latter 
quantity. This effectively turns a problem of definite integration into one of indefinite integration 
(in a different set of variables) . We shall use this approach to give concrete expressions for all the 
box integrals. It is also possible, as we shall see in the next section, to derive a purely algebraic set 
of equations for the n-point integrals /„, in which a new unknown quantity enters only at 0(e). 



3. Algebraic Equations for n-Point One-Loop Integrals 

In this section we will derive a set of algebraic equations for the general w-point one-loop 
integrals. Some of the equations are of use in the partial differential equation approach of section 4; 
others can be used to determine the n-point scalar integrals for n > 5 in terms of box integrals, 
in an entirely algebraic fashion (subject to some subtleties for n > 7, which are explained in 
appendix VI). The equations have been derived in ref. [10] using a momentum-space representation 
of the loop integrals. Here we will derive the same general equations using the Feynman parameter 
representation; in this derivation the equations arise from the consideration of integrals of total 
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derivatives of the Feynman parameters.* For a specific, simple example of the following general 
derivation, we refer the reader to the beginning of section 4. 
The total derivatives we will consider are 



Jn;m = r(n - 3 + e) / da n -i 

Jo Jo 

d 1 



1 — a„_i pi — oi— <Z2 a m ctn-i 

da n - 2 ■ ■ ■ / (ia r , 
■/ o 



da r 



^2i,j = l Sij a i a j 



n — 3+e 



(3-1) 



a„ = l — a!— a 2 a„_i 



There are two ways to evaluate J n;m . First, one can carry out the differentiation with respect to 
a m , to get 



J n 



-2T(n-2 + e) / cfX <f(l - 



S 7" = 1 ( "^"l j a j ^nj a j ) 



Si, j = l Sij a i a j 



n-2+e 



(3.2) 



= -2 j JJ a € ) ^ I n [(S mj - S nj ) a j aj] . 



l=i 



Second, one can perform the integral over a m . At the lower integration endpoint, a m is set to 0. 
The remaining (n — l)-point integral corresponds to removing the propagator parametrized by a m 
- i.e., the propagator between lines (m — 1) and m — from the original n-point (scalar) integral; 
we denote such a "daughter" integral of I n (J n ) by Similarly, at the upper integration 

endpoint a n is set to 0, yielding the (n — l)-point integral I^-i- It is always possible to choose the 
cti variables for the integrals so that they are the same as those for the parent integral I n . 
Having made this choice, the second evaluation of J n;m gives 

= /&[!]- Ji-iM = (flat 



Jn 



\l=l 



j(») 



p(m) 



n—1 n—1 



a r 



(3.3) 



Equating (3.2) and (3.3), using Sij = pij/(aiOtj) and the definitions (2.15) of the reduced integrals, 
and relabelling the index m — > i, we have 



E 

i=i 



Pij Pnj 
Oii C^rj. 



4 [«j 



f{i) n 

1 n-l 1 n-l 



tin) 



OL r . 



i = 1, 2, . . . , n — 1. 



(3.4) 



* The motivation for considering such objects arose from the observation that the field-theory limit of integrals 
of total derivatives in string theory yields expressions that are sums of loop integrals with differing numbers of 
external legs (multiplied by various coefficients); these sums must necessarily vanish because the world-sheets in 
the string loop expansion have no boundaries, when appropriate analytic continuations of the external momenta 
are used. 

t For more explicit examples of this notation, see the beginning of subsection 4.2. 
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We would like to solve for the n one-parameter integrals I n [dj]- To do so we supplement the n — 1 
equations (3.4) with the equation that follows from the constraint on the Feynman parameters, 
Sj=i a 3 = !> namely 

n 

^OLj i n [aj\ = I n [l] = I n - (3.5) 
i=i 

Before solving equations (3.4), (3.5), we introduce a little more notation and some "kinematic" 
results from ref. [10]. We define the Gram determinant of the (n — l)-vector system associated with 
the n-point integral by 

A n = det'(2jfei • kj), (3.6) 

where the prime signifies that one of the n vectors ki is to be omitted before taking the determinant; 
due to momentum conservation, ^ ki = 0, any one of the vectors may be omitted.* Next we 
introduce the rescaled Gram determinant, 

n 

A„ = (n«<)An, (3.7) 

e=i 

which has a simple bilinear representation in terms of the variables a^. 

n 

A„ = ^2 Vij^iaj. (3.8) 

Here r/ij is independent of the a^; in fact r\ is proportional [10] to the inverse of the matrix p defined 
in equation (2.10): 

p = N nV -\ 7] = N n p~\ N n = 2 re - 1 detp. (3.9) 

We also define the variables 7^ by 



= ^ i dA n 

i=i 

They are in a sense conjugate to the a>i variables: 



(3.10) 



pij fixed 



If we define 



^PijKj = N nOti- (3.11) 
i=i 



Rki = Vki ~ ^p- , (3.12) 



$ The notation for, and normalization of, the Gram determinant in equation (3.6) differ from other conventions 
in the literature, e.g. references [5,18]. 
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then we may note the following identity, 

n 

Y^UiRki = 0, (3.13) 

i=l 

which follows from equations (3.8) and (3.10). Thus if we multiply both sides of equation (3.4) by 
otiRki, and then sum over i, the terms not containing a, will drop out, leaving us with 

n / \ 1 n 

[VkiPij ~ j^HPij) L[aj] = TjJ^RkJnli , (3-14) 
.,7=i V A n / ^ . =1 



or 



A^nM = \ Y. R^n-l + ^ n £ . (3.15) 

i=l ^ n j = l 

Performing the sum on the right-hand side with the help of (3.5), dividing by N n , and writing 
out the definition of Rki, we obtain 

^K] = ^E(*-^)a + fjn. (3.16) 

Combining this set of equations with the derivative representation (2.20) for m = 1, we obtain a 
system of partial differential equations for the n-point scalar integral, 

1 dL 



^E(*-F)^ + P"- P.!7) 



n - 4 + 2e 2iV„ ^ V 111 A n > n ~ L A n 

Section 4 is devoted to solving these equations for various scalar box integrals. 

In ref. [10] a momentum-space representation was used to derive an algebraic equation that 
involved only scalar integrals, at the expense of introducing a new object, In =6 ~ 2e - The object 
I® =6 ~ 2e comes from an integral in D = 4 — 2e with two loop-momenta inserted in the numerator, 
but it can also be interpreted as the n-point scalar integral in two higher dimensions. The latter 
interpretation is helpful for understanding the properties of I® =e ~ 2e as e — > 0, which are needed in 
order to use the "dimension-changing" equation to obtain D = 4 — 2e scalar integrals through 0{e ). 
We shall now re-derive this equation using Feynman parameter representations of the integrals. 

The integral I^ =6 ~ 2e [l] is most easily obtained from the D = 4 — 2e equation (2.2) by letting 
e — ► e — 1, 

/f= 6 - 2e [l] = T(n-?, + e) f\ n ai 5{l-Y.^)— —3+7, (3.18) 

Yli,j = l Sij a i a j 

It may also be obtained by inserting one power of the scalar denominator of the D = 4 — 2e integral 
into the numerator (summations are implicit in the following derivation): 

^ =6 - 2 11] = r(n-2 + c) 7 " =4 " 2£[ ^ Q ^ ] = ^JT~e Plj =4 " 2e WoiH'iM " ( 3 ' 19 ) 
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In terms of the reduced integrals (2.15), and using the derivative representation (2.20), we have 



fD=6-2e 



(n - 3 + e) (n - 4 + 2e) (n - 5 + 2e) Pij dc^da 



(3.20) 



We may now evaluate the right-hand-side of (3.20) using equation (3.17) to replace the deriva- 
tives, and also using the relations between p, r\ and 7 to simplify the expressions: 



n-A + 2e Pl3 'da l da 



d 2 h 



Pij 



1 



da f \ 27V, 



»7jfc 



7j7fc 



f(fc) ■ Is f 
A~ 



2JV, 



Al 



A, 



n-l + 



1 



2N„ 



Vjk 



A n 



dl 



Pij- 



(k) 

n-l 



da. 



+ 



pij 

A„, 



Vij ~ 2 



A n 



n-l Ik f(k) 1 
2 A„ - 1 + 2 



I n + (n - 4 + 2e) 
_<9 7fc g 





7i7fc 


Giv n 


Vik x 

A n 



-'n-l + 



7i 



e ) i+((^-2) + (n-4+ 2e ))^/ 1 



Now I^-L is actually independent of (since cifc has been set to in I n K -i)', also 

^ft-'i = (n-5 + 2c)e ) i- 



da,- 



(3.21) 



So we obtain 



n - 4 + 2e Plj 'daidaj 
which can be solved for /„ using equation (3.20), 

E7.£i + (n-5 + 2 e )A n # 



= (n-3 + e) 



2AL 



_7fc_#(k) _ 
fc=i ^™ ^ 



(3.22) 



1 

2iV7 



=6-2e 



i=l 



(3.23) 



In ref. [10] it is shown how to use this equation to obtain n-point integrals with n > 6. However, 
for n > 7 there are some complications, which are discussed in appendix VI. In this paper our main 
interest is the pentagon integral (n = 5). For the scalar pentagon integral it suffices to note that 
the integral j^= 6 ~ 2e [ s finite as e —> 0, because the D = 6 scalar pentagon integral possesses neither 
ultraviolet nor infrared divergences (soft or collinear), and also that the coefficient of jD= 6 ~ 2e [ n 
equation (3.23) is of order e. Therefore to 0(e ) the general scalar pentagon integral is given by 
the sum of five scalar box integrals, 



27V, 



-J>/f + O(e) 

^ i=i 



(3.24) 
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A schematic depiction of this equation, with the coefficients suppressed, is given in fig. 1. For the 
tensor pentagon integrals we have to keep the f^= 6 - 2e term around a while longer (see section 6). 

One further equation for all n can be obtained by eliminating I n from equation (3.16) using 
equation (3.23), with the result 



1 



dl„ 



n-4 + 2e <9a; 



= I„.\aA = 



2N n 



E^*-i + (--5 + 26) 7l /^ 6 - 2e 



(3.25) 



Since the D = 6 scalar box is also finite, setting n = 4 in equation (3.25) yields a simple set of 
partial differential equations for the box integrals, through 0(e): 



dh 

don 



e 

iV4 



^Vijl? + (-l + 2e) 7l /4 



D=6-2e 



^E%^ + O(e). 



(3.26) 



i=i 



The right-hand-side depends only on the infrared singular pieces of the triangle integrals. 

This completes our re-derivation of general all-n results presented in ref. [10]; we now apply 
these results to various box and pentagon integrals. 



4. Partial Differential Equation Technique 

In this section, we solve the partial differential equations (3.17), (3.26) for scalar box integrals 
with all internal lines massless, but with 0, 1, 2 or 3 massive external lines. 

4.1 The Massless Box Integral 

We begin with the box integral with all external lines massless, 

h[l] = T(2 + e) / d 4 a, «5(1-E^)t \ ^ ■ t 4 - 1 ) 

Jo [— said3 — ta2di\ 

This integral is simple enough to perform directly after the following change of variables [19] which 
factorizes the integrand^ 

a-i = y0-~x), a 2 = z(l-y), a 3 = (1 - y)(l - z), a 4 = xy. (4.2) 

However, our purpose here is to illustrate the partial differential equation technique, including the 
derivation of the equations, via this simple example. 

' J. Vermaseren has pointed out to us that the factorization of the integrand in terms of x, y, z arises naturally 
if one combines pairs of propagators using Feynman parameters, and then combines the two resulting factors 
using another Feynman parameter. See also ref. [19]. 
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As noted above, algebraic equations for Feynman parameter integrals can be obtained by- 
considering integrals of total derivatives. Here we consider the box integral h[l], with the parameter 
04 eliminated, and with the integrand differentiated with respect to a\\ 

pi pl—a 3 pi— a 2 — a 3 q ^ 

J 4;1 = r(l + e) / da 3 / da 2 / da, — - — . (4.3) 

Jo Jo Jo oa i [— sa,a 3 — ta 2 (l — 0,1 — 0,2 — a 3 )\ 

Observe that J A; i can be evaluated in two ways, either by explicit differentiation, or by evaluating 
the integrand at the boundaries 04 = 1 — a, — a 2 — a 3 = and ai = 0. The boundary terms yield 

r(l + e) f da 1 da 2 da 3 5{ l l ~^= 1 ^ ) - T(l + e) [\a 2 da 3 da j^ ~ , (4.4) 

Jo [-saia 3 J 1+e J [-ta 2 a A \ 1+t 

which is the difference of two triangle integrals, each with one massive external leg, as depicted in 
fig. 2. These integrals are easily evaluated, 

Il-{s) ^ F(l + e) fd*a ^-^=^ = , (4.5) 

J [sa 1 a 3 \ 1 + e e 2 

where 

r(i + 6)r 2 (i- £ ) 

rr ~ T(l-2e) 

is a ubiquitous prefactor. Thus 

J 4;i = ^((-s)- 1 - 6 -^)- 1 - 6 ) • (4-6) 
The other way of evaluating equation (4.3), explicit differentiation, yields 

J 4; i = -r(2 + e) r^a^l-E^) " Sa3 + ta2 
Jo 



[— saia 3 — ta 2 a4] 

= h[sa 3 -ta 2 ] = ^\Y\ a i 



2+e 



i=l 



Id Id 
+ 



a, da 3 CX4 da 2 



(4.7) 



h 



where we have used equations (2.11), (2.15) and (2.17) in the last step. 

Equations (4.6) and (4.7) together constitute one differential equation for I4. In fact, due to 
the symmetries of the original integral, total derivatives in other Feynman parameters do not yield 
independent equations. Instead, we recognize at this stage that I4 is really a function of s and t 
alone, not of all four c^, 

h = h(s,t) = / 4 (-(a 1 a 3 )" 1 ,-(a2a4)" 1 ), (4.8) 

so that 

1 dI A _ 2 dh J_£^£ _ + 2<H± /a q\ 



«i da 3 ds ' 04 da 2 dt 
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Combining equations (4.6), (4.7) and (4.9), we see that l4(s,t) satisfies the partial differential 
equation 

s2 f - < 2 § - -¥ st [ ( - s) " 1_! - ( - t)_1 " ! ] ■ (4 - io> 

We still need one additional equation, which comes from the fact that the dimension of I4 is equal 
to — e x dimension ( s, t), so that 

8 ar + = - £/ * • (4 - u) 

Equations (4.10) and (4.11) form a complete set of partial differential equations. 
If we consider instead of I4 the dimensionless quantity I®, defined by 

/ 4 °M)=(-^) h(s,t), (4.12) 
we see that it is a function only of the ratio x = t/s, an d that 

In terms of x, the first equation (4.10) becomes 

dil _ 2r r (x'-X' 1 ) 



(4.14) 



rf X e (l + x) 1+e 

One can solve this differential equation to all orders in e as follows. We observe that the 
transformation x X -1 interchanges the two terms on the right-hand side. Taking the second 
term, shifting x ~~ * X ~ lj using the hyper geometric function formulas 

f z c+1 
dz z c p F q ({ai};{bi};z) = — — p+ iF q+1 ({at}, c + 1; {&;}, c + 2; z) (4.15) 

and 

iF (Z;z) = (l-z)S, (4.16) 

the hypergeometric function identity 

2J F 1 (l,-e;l-e;l + X ) = (~xY 2 i ? i(-e, -e; 1 - e; 1 + x _1 ) , (4.17) 

and using the interchange of x an d X -1 to furnish the first term, we obtain (note that x should be 
thought of as having a small imaginary part in order to avoid difficulties with branch cuts) 

fo _ 2rr 



(1 + x" 1 ) 2^1 (1, -e; 1 - e; 1 + x" 1 ) + (1 + x)" e 2^1 (1, -e; 1 - e; 1 + x) 



^ (1 + x" 1 ) 6 [(-X" 1 ) 6 2^1 (-C -6; 1 - e; 1 + x) + (-x) £ X _£ 2^1 (-€, -e; 1 - e; 1 + x" 1 )] 

(4.18) 
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The constant of integration may be determined by evaluating the integral (4.1) directly at some 
convenient value of x, say x = t/s = 1 (s = t = —1), 



I 4 °(X = 1) = 2- e / 4 [l] = 2"T(2 + e 



1 4 <5(1-Eai 



d o,- 



[aia 3 + a 2 a 4 ] 2+e 

= 2"T(2 + e) [ dx [ dy [ dz y'^^l - y)' 1 '" [(1 - x)(l - z) + 
Jo Jo ■/ 



(4.19) 

where we have made the change of variables (4.2). The y and z integrals are elementary and leave 
us with a standard hyper geometric integral, 



-2-e 



#(X = 1) = - 



2 1 - £ r r 



/' 

Jo 



dx 



x 



-l-e 



(1-x)- 



-l-e 



e Jo 1 — 2x 

2 1 " £ r r r(-e)r(l) 



e T(l-e) 5^o V 



lim ( 2 Fi(l, -e; 1 - e; 2 + z<5) - 2 Fi(l, 1; 1 - e; 2 + z<5) ) (4.20) 



2 i " e r r 



lim 

5^0 



( 2 iq(l,-e; l-e; 2 + iS) + 2 F 1 (l,-e; l-e; 2 - i<5)) . 



Comparing with the first line of equation (4.18), we see that the constant of integration vanishes. 

Alternatively, we may solve equation (4.14) order by order in e. Observe that / 4 must con- 
tain 1/e 2 poles from the overlap of collinear and soft singularities. As the right-hand side of the 
differential equation only contains a single power of 1/e, this leading pole should be multiplied by 
something to the ±e power, so that one power of e is cancelled upon differentiation. Through O(e ), 
we then have 



II = r r ||[(l + x)- £ + (l + X- 1 )- £ ] -\n 2 x-n 2 } + 0(e 
(— ) H 



(4.21) 



" + O(e) , 



where the constant of integration can be fixed as in the all-orders solution. 

Restoring the prefactor (— (s + t)/st) £ , and expressing the result in terms of the aj, we have 



h[l] 



2r r 



a 2 «4) e 2 Fi I — e, — e; 1 — e; 1 



aids 
0.201a 



+ {-oliolzY 2 F\ -e,-e; l-e; 1 



a 2 a 4 \ 
evicts) 



— ((aia 3 ) e + (a 2 a4) e ) - In 



2 / aia 3 \ 



7T 



+ 0(e) 



(4.22) 

In this form, the differentiation formula (2.17) may be applied to the scalar integral I4 to obtain 
the integrals with arbitrary Feynman parameter polynomials inserted. Because of the T(l — m + 2e) 
prefactor in (2.17), the 0(e) terms in I4 contribute to the polynomial integrals at O(e ). Instead 
of displaying the 0(e) terms in I 4 explicitly, we quote the reduced integrals with one parameter 
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inserted, Z^Oj], to O(e ): 



aih[ai] = a 3 / 4 [a3] = r r <j \(a 2 a 4 ) e - - 



0103 



2 V «i«3 + 0204 



In 2 



204/ J J 



a 

020:4 



+ 0(e), 



02/4(02] = 04^4(04] = r T \ — (aia 3 ) f 



0204 \ 



2 \a1a3 + 0204 



+ 0(e). 

J [ \a2a4J J J 

(4.23) 

The latter integrals may be differentiated further to obtain through O(e ) the integral with any 
polynomial of the Feynman parameters inserted. 

As mentioned previously, the branch cuts can be obtained by inserting the ie associated with 
each kinematic variable, 



(4.24) 



(-s)~ e |sr e e + ^ e0(s) , 
ln(— s) — > ln|s| — iir@(s) , 
where @(x) is the usual Heavyside function: Q(x) = 1 for x > and 0(x) = for x < 0. For the 
massless scalar box we therefore obtain 



I 4 [1] 



rr 1 



(47r) 2 st [ e : 



V 47T/U 2 / 



-In 2 



+ 2m(@(s) - 9(t))ln 



7T 



i_ ( e (s) -G(t)) 2 +0(6), 



(4.25) 



where s and i are the Mandelstam variables defined below equation (2.6). 

4.2 The Box Integral with One External Mass 

Following the same techniques, we can obtain partial differential equations for boxes with one 
external massive leg (or equivalently, one external leg off-shell), 

1 



I 4 m (s 1 ,s 2 ,m 



T(2 + e 



f d±a 5(1 - Y,i 
Jo 



n 2+e • 



(4.26) 



[— saias — ta 2 ai — m\a4,a\] 

(This integral could also be evaluated using the same change of variables (4.2) as for the massless 
box.) By analogy with equation (4.4), such integrals will clearly arise in the consideration of 
massless pentagon integrals. Following the conventions of section 3, we label these boxes by 1^ 
when the momentum invariant Sj_i ; j for the adjacent legs (i — 1) and i of the pentagon diagram 
serves as the "mass" of the massive leg of the box. For example, 

4 5) [1] = / 4 lm (si2,S23,S4 5 ) 

1 (4.27) 



= r(2 + e) f d 4 a5(l-Z t ^) — 
Jo \—i 



[-SI2O1CI3 - S230204 - S45040l] 2+e 

is the box integral arising from the diagram depicted in fig. 3, in which a tree with external legs 4 
and 5 is attached to a four-point loop. Note that the scalar denominator for the integral (4.27) 
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can be obtained from the massless pentagon denominator by setting the parameter 05 to zero. 
Similarly, 1^ can be obtained by setting a« — ► in the massless pentagon. 

From these remarks it is clear that the change of integration variables described earlier for 
the pentagon can be used here to remove the kinematic factors from the denominator of the box 
integral, 

4 5) [1] = r(2 + e)(n°* 



•L 4 <f (l-£«i) (E^iOjtij) 
,- d U 7TT- 

'0 [uiu 3 + u 2 u 4 + U4U1] 



2e 



(4.28) 



where a, are given by equation (2.13). (These variables a, should not be confused with the corre- 
sponding a, for the massless box.) The other integrals that will arise, 



r(i) 



- ^4 m ( S 23> s 34, S 51 ), < , 



/i 2 ' 1 — / 4 m (s3 4 , S 4 5, S12), 



r(3) flm/ \ 
i 4 — i 4 (S45, S51, S23J 



4 4) = il m (351,812,334), 



(4.29) 



(5) "CO 

can be obtained from I\ by cyclic permutation of the CKj. We define the reduced integral, I\ or 
/ lm , via equation (2.15). 

We can now apply the general results of section 3 to the example of equation (4.27). The 
matrix p defined in equation (2.10) is now given by 



,1m 



/0 1 1\ 

1 

10 

Vi 1 0/ 



(4.30) 



so that N\ m = |, and the rescaled Gram determinant is given, using eqs. (3.8) and (3.9), by 



A4 111 = 2(«ia 3 + a 2 a 4 - a 2 a 3 ) 



(4.31) 



Using equations (3.8) and (3.10), the explicit values of the quantities 7/™ and rjj™ can be read off 
from (4.31): 



7x lm 



«3, 



«4 - «3, 73 m 



,1m 



n lm 
742 



Oil ~ OL2, 74 



«2, 



remaining 77, 



lm 

ij 



0. 



(4.32) 



In terms of these quantities, the differential equations (3.16) read 



dl 



lm 



dctj 



= 2e 



= 2e 



lm lm\ lm 



V Tl lm - ' ' J f j) + U J 1 



i=i 
4 



A 



! ^A- 

■ 13 i , - J/ 



.7 = 1 



4 doiidctj 3 



<9a 7 - 



(4.33) 
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The triangle integrals appearing on the right-hand-side of (4.33) include both the triangle integral 
with one external massive leg, Ig m , defined in equation (4.5), and the triangle with two external 
masses, 



r r (si) e - (-s 2 )' 



6 2 (Si) ~ (-8 2 ) ■ 

Explicitly, the following reduced triangle integrals appear (defined again via (2.15)): 



j(l) _ 


1 


rim 
J 3 




f(2) 


1 


r2m 
J 3 




f(3) . 


1 


r2m 
J 3 


a\OL 2 OL^ 


j(4) _ 
J 3 — 


1 


rim 
J 3 


aia 2 a 3 



1 jlm ( - 1 \ = r r «4«2 

*3"4 3 V a 2a4/ e 2 «3 



ai«3 «4Qi 



ai« 3 



e 2 a 3 


5 






rr 




)- 


c 2 


a\ 




rr 




)- 


c 2 


«4 








e 2 a 2 







CK4 — «3 



«i — a 2 



(4.34) 



(4.35) 



The differential equations (4.33) have the solution, 



fim _ 2r r 

J 4 - — o- 



. '»;(!.•, • n 2 n i n-jn.A 

(-a 3 (ai-a 2 )) 2-fi -e, -e; 1 - e; ? r 

« 3 (qi - a 2 ) / 



+ (-a 2 (a 4 - a 3 )) e 2 Fi ( -e, -e; 1 - e; 



aia 3 + a 2 a 4 — a 2 a 3 \ 
a 2 (a 4 -a 3 ) / 

e / aia 3 + a 2 a 4 - a 2 a 3 \ 

- ((ai - a 2 )(a 4 - a 3 )j 2 iq -e, -e; 1 - e; -— — 

V (ai - a 2 )(a 4 - a 3 ) J 



2r r 



(Aim 



(A lm 
-e, -e, 1 - e, j^if^ 

-(7W^(- £ ,- £ ;i-e;-^ 
= 2rr (^ + Li2 f 1 _^) +Li2 (' 1 _2l^ 

a 2 / V «3 



+ 0(e) 



(4.36) 



where Li 2 is the dilogarithm [20], which satisfies 

d 



dx 



Li 2 (l-x) = 



ln(x) 

1^ ' 



(4.37) 



and also the identity 



Li 2 (l-x)+Li 2 (l-x- 1 ) 



= -2 ln2 ( x )' 
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x > 0. 



(4.38) 



In principle, one could also add any solution of the homogeneous equations, (4.33) with 1^ set to 
zero. The coefficient of such a solution vanishes, as may again be demonstrated by evaluation at a 
special kinematic point. 

In terms of momentum invariants, the unreduced integral is 



r( 5 ) 



(-S12) e (-s 2 3)~ 



•S12S23 



(-S45)" 



+ Li. 



i 2 (l-^)+Li 2 (l- 



£23 \ _ 
S 45 / 6 



+ 0(e), (4.39) 



or alternatively, after using the dilogarithm identity (4.38) and rearranging the terms 



r(5) _ 
J 4 — 



S12S23 
2 Li 2 



(-Sl 2 )" £ + (-S2 3 ) 



-S45J 



1 



^ ] - 2 Li 2 

S12 



S45 



■523 



«12 



1 _ _^ _ i n - 1±± _ _ + (e). 



$23 



(4.40) 



This second form is appropriate for studying the limit S45 — ► as we do at the end of this section. 

Including the overall normalization factors appropriate for the momentum-space integral (2.1) 
yields 

2 .(p|)-7-£2|)-'(-fl|) , [I + Li 2 ( 1 -fli) +Li2 ( 1 _f2i)4' 
V 47r u / \4iru 2 J V Anur J e 2 V s 45 / V s 45 / 6 



r(5) 



(47r) 2 si 2 s 23 \47r^ 2 



+ 0(e), 
(4.41) 

in agreement with the results of refs. [21]. The correct analytic continuation to the physical region 



+ ie. 



can be obtained from this expression by taking s 

As with the massless box, it is useful to quote the integrals / 4 m [ai] to 0(e ); further differ- 
entiation of them will give any desired integral to 0(e ) as well. The I 4 m [ai] may be read off 
from equations (4.33) and (4.36). We rewrite them in terms of a combination of dilogarithms and 
logarithms that will reappear in the massless pentagon tensor integrals: 



4 5 Vi] = 

4 5) N = 



r T 



1 a\{a\ - a\) 



a 2 



0L2 



+ 



013 L 5 



aias + a 2 a 4 — a 2 «3 



+ 



+ 



a%{a\ 



cci — a 2 
at\{a\ - a £ 3 ] 



where 



Li = Li? 



r r 



03 a 4 
1 a\{a% - a\) 



«3 



(a 4 - a 3 ) L 5 



ot\oi3 + a 2 a 4 — a 2 a3 
(ai - a 2 ) L 5 



«i+2 



a 4 — «3 



+ Li 2 1 



+ 



ai« 3 + a 2 a 4 
"2 £5 



a 2 a 3 



(4.42) 



«ia 3 + a 2 a 4 — a 2 a 3 



Ui-l 



a, 



+ In 



Oti+2 



In 



«i-2 



7T 

~6 



(4.43) 



Note that L, vanishes as ai+iai^ + a i+2 Qj_i — ai +2 oti- 2 — > 0, so the /^[a,] are not singular in 
that limit. 
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4.3 Box Integrals with Two External Masses 

In order to evaluate the pentagon integral with one external mass, or the all-massless hexagon 
integral, one needs box integrals with two external masses, of which there are two types, which 
we will call 'easy' and 'hard'. Both of these integrals have been performed previously [22,9]. The 
'easy' box, with external masses at diagonally opposite corners, can be done with the same change 
of variables (4.2) described in section 4. We will not discuss it further, but merely quote the result, 



f 2me [l] = T(2 + e 



d A ai 



6{1 -£"0 



,2+e 



2r r 



st — m 



Li 2 (l-^ 

s 



[— saia3 — ta 2 a 4 — m\a\a2 — m\a^a^ 

R {i «-)- + <-«r - (-»?)- - <--»-•) + u, (i - 



Lb 1 



m 



Lio 1 



m\ 



Lb 1 



^)-^ 2 (f)} 



The 'hard' box, with external masses at adjacent corners (legs 3 and 4), 



j2mh 



[1] 



T(2 + e)^ d 4 a- 



8(1 ~ 



n 2+e 



+ 0(e) 
(4.44) 

(4.45) 



[— sa,\a,3 — ta 2 a 4 — m^a^a^ — m\a 4 a\\ 
cannot be easily done this way; but it is amenable to the partial differential equation technique. 
We change to «j variables defined by 

1 



1 



s = — - 



t = -- 



aia 3 

Then the matrix p is given by 



a 2 a 4 



m-i = —- 



m 2 



1 



a 4 a\ 



„2mfi 




(4.46) 



we have iV| m/l = ^ and the rescaled Gram determinant is 



A| m = 2(ai«3 + ct 2 a 4 — aiQ2 — a 2 a 3 + a 2 ) 



(4.47) 



from which 7 2mh and J? 2 ™' 1 can be obtained via equations (3.8) and (3.10). 

To obtain the box integral to O(e ), we can use the simple partial differential equations (3.26), 
which are sensitive only to the pieces of the triangle integrals that are singular as e — > 0. In 
particular, the three-mass triangle does not contribute, because it is finite as e — ► 0. We find that 
to O(e ), 



d f2mh 



da.; 



J2iah f{j) 



2r r 



1 r) 2mh 1 

e a 2 



2mh ln ( a 2/a 3 ) 



a 2 — «3 



„2mft ln ("i/ a 2) , ^ 2mfa ln(ai« 3 ) 

%3 ~~ " r T] u 



(4.48) 



«1 — «2 



a 2 
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Writing 



f2m/i 

4 



4r + -X-! +X + c 



(4.49) 



and solving the differential equations for X-\ and Xq, we find 
X_i = 21na 2 , 



or 



«2 



2Li 2 1- — +2Li 2 1-— + 21n"a 2 , 



«3 



a 2 



1 



+ 2Li 2 1 



«2 / 



+ 2Li 2 1 



«3 

a 2 



+ c 



+ O(e) 



(4.50) 



(4.51) 



The constant Co may be determined by computing the function at a specific point, say where 
all the cc, are equal; the resulting integral is evaluated explicitly in appendix IV, whence we find 
Co = 0. Finally, rewriting the result (4.51) back in terms of the conventional kinematic variables 
yields 



2\e 



ir h [l]=r r 



l+e 



1 



+ 2Li 2 1- 



+ 2Li 2 1- 



+ 0(e). (4.52) 



Using the dilogarithm identity (4.38) and rearranging the terms this can be written in the 
alternative form 

1 {-ml)- e {-ml)- e 



j2mh 
1 A 



[i] 



rr 
st I e 



,2\-e 



(-«)-« + (-*)-«- (-m^)-' - (-mfl 



+ 



(-«)" 



2Li2 ( 1 _^)_ 2Li2 ( 1 _^)_ ln2(?) } +0(e) , 



(4.53) 



which is more convenient for studying the massless limit, as we do at the end of this section 

4.4 The Box Integral with Three External Masses 

Here we compute the three-mass scalar box integral, 

<y(i-E«i) 



/ 4 3m [l] = r(2 + e) J d 4 



[— said3 — ta 2 a4 — m^a^a-i — m^a^a^ — m\a^a\\ 2+6 
We again use the partial differential equations (3.26), with the change-of- variables 



(4.54) 



s = — - 



t = -- 



A 



a 2 «4 



m = — - 



mo = — - 



m„ = — - 



a 2 a 3 



(4.55) 



which is the same as that used for the hard two- mass box, except that now A / 0. The matrix p 
becomes 



3m 

Pa 



1 1\ 
1 I A 1 

1 A 1 
1 1 1 oJ 



(4.56) 
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the normalization factor in (3.26) is N^ r 

2A 



= 5(1 — A) 2 , and the matrix 77 used to construct Af 11 is 



V 



3m 




-l-A 
2 

-(l-A) 
1 - A 



1 - A 
-(l-A) 







-A(l-A)\ 
l-A 


/ 



(4.57) 



We expand if™ as 



f 3m 
J 4 



ci(A) 



+ X (a t ,X) + co(A) +0(e) 



(4.58) 



(We will see below that there is no 1/e 2 singularity.) To solve the partial differential equations 
order-by-order in e we need to first know ci(A). We know that ci(A) is independent of the 
because the daughter triangles here are of the two-mass and three-mass varieties; the two-mass 
triangle has a 1/e pole, which feeds into Xo(a>i,\), while the three- mass triangle is finite and can 
be ignored altogether. So we may compute ci(A) by doing the integral lf m for the special choice 
of all ai = 1. We should compute the finite part of the integral while we're at it, since this result 
will fix the constant of integration c (A). This computation is done in appendix IV, where we find 

In A 



Cl (A) 



A 



(4.59) 



Next we solve the partial differential equations (3.26). Plug the expansion of tf m (equa- 
tion (4.58)) and the divergent pieces of the 2-mass triangles 1^ and I^ 4 \ 

r r ln(ai/a 2 ) 



?(3) 



r(4) 



e a>i — a 2 
rr ln(Aai/a2) 



+ O(e ) , 



(4.60) 



+ O(e ) , 



e Aai — C12 

into the far right-hand-side of (3.26) and use the result (4.57) for 77^, to get 



dX 
da 1 
dX 
da 2 
dX 
da 3 
dX 
804 



1 - A 
2 

1 - A 
0, 

. 



\n(ai/a 2 ) + ^ln(Aai/a2) 
a i — (*2 A«i — «2 

ln(ai/a 2 ) ln(Aai/a 2 ) 
a\ — «2 A«i — «2 



(4.61) 



Solving these equations for X (ai, A), and fixing the constant Cq(A) using equation (IV.6), yields 



f3m 

J 4 — 



A 



In A 



+ 2 Li 2 1 - 



a 2 



- 2 Li 2 1 - 



A«i 

a 2 



1 



+ 0(e) . 



+ 2 Li 2 (l-A) + 21nA lna 2 - - In A 



(4.62) 
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Returning to the original kinematic variables, and using dilogarithm identities [20], we get 

^M,™?) = . Z s- 2 { | [(~*r + - (-mi)- 5 - (-mg)-' - (-ml)-" 

1 (-m|)- £ (-mj)- £ , 1 (-m|)- £ (-ml)- e 
e 2 (-i)" e e 2 (-s)" e 

(4.63) 

4.5 The Box Integral with Four External Masses 

The four- mass box integral is infrared finite and has been performed in D = 4 in ref. [6]; a 
compact expression is given in ref. [15]. Amusingly, the partial differential equations (3.26) for it 
are trivial, because the three-mass triangles appearing on the right-hand-side are non-singular as 
e — > 0. In other words, through 0(e), the reduced four-mass box cannot depend on the a,, but 
only on the two other, dimensionless variables, say Ai and A 2 , where we define 

I i 2 Ai 2 A 2 2 1 2 1 

m 1 = , m 2 = , m 3 = , m 4 



a\a 3 a 2 a^ ol\ol 2 " a 2 a 3 a 3 a^ a^a\ 

(4.64) 

One can check that the answer Do(s,t,m?) given in ref. [15] does have this property — when the 
integral is divided by a\a 2 a 3 a^ it depends only on Ai and A 2 - Indeed, 

It = = I/ L i 2 (i(l-A 1 + A 2 + r)) - Li 2 (|(l-A 1 + A 2 -r)) 

1 1 a i r I 

+ Li 2 (2^(1 - Ai — A 2 — r)) - Li 2 (^(1 - Ai — A 2 + r)) (4.65) 



where 



r = y/l - 2Ai - 2A 2 + Xj - 2A X A 2 + A| . (4.66) 



4.6 The Massless Limit of Massive Boxes 

In general there is no reason for the massless limits to be smooth. The limit of taking a mass 
to zero does not necessarily commute with the 1/e expansion of dimensional regularization, which 
has been truncated at O(e ). For e < (as is required to regulate the infrared divergences in the 
box integrals), we see that the single external mass box ll m (s, t, m|) (given in equation (4.40) with 
s = S12, t = s 2 3, m\ = S45) goes over smoothly to the massless box l2 m (s,t) as m 4 — > 0, and 
the easy two-mass box /| me (s, t, m§, m|) goes over smoothly to the one-mass box /| m (s,t,m 2 ) as 
m 3 -► 0. On the other hand, the limits, /| m/l -► 7 4 lm , /| m -► /| me , and /f m -► /| m?l , are not 
smooth: in each of these cases there are "missing" dilogarithms. 
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The fact that some of the above limits happen to be smooth, with only the exponentiation 
of the logarithms (— s)~ £ , (—t)~ e , (—mf)~ £ , can be understood from the representation (3.23) (for 
n = 4) of the D = 4 — 2e box integral as the sum of D = 4 — 2e triangles and a D = 6 — 2e box 
integral. The D = 6 — 2e box integral is infrared (and ultraviolet) convergent for any choice of 
mass, so it has a smooth limit as any mass goes to zero. The D = 4 — 2e triangles appearing in the 
representations (3.23) for I®" 1 , l\ m and l| me have either one or two nonvanishing external masses; 
these integrals can be written in closed form to all order in e merely by exponentiating logarithms. 
(See equations (4.5) and (4.34).) In contrast, the representations (3.23) of the box integrals /f m/l 
and /| m require the triangle with three external masses, /| m , whose all-orders-in-e form (V.ll) is 
considerably more complicated, involving hyper geometric functions. One should not expect that 
these latter box integrals, truncated to 0(e°), could be made to have smooth limits simply by 
exponentiating logarithms. 



5. Algebraic Approach to Pentagon Integrals 

It is possible to solve the partial differential equations (3.17) for the massless scalar pentagon 
through O(e ). However, a simpler approach, which works equally well for arbitrary pentagon 
kinematics, is to use the general algebraic equation (3.23) derived in section 3 to express the scalar 
pentagon integral I5 as a sum of five scalar box integrals, up to 0(e) corrections: 



1 

2N~ 5 



(0 



+ 2eA 5 / 5 D = 6 - 



2t 



i=l 



(5.1) 



(See also fig. 1.) To give an explicit expression for the pentagon, we need only collect the relevant 
scalar box integrals from section 4, and compute the kinematic coefficients N$, A5, 7$ and rjij. 
(The rjij are relevant for computing tensor integrals.) We now do this for the all-massless pentagon 
integral, and for the pentagon with one external mass. 

5.1 The Massless Pentagon Integral 

For the massless pentagon, equation (2.7) for the scalar denominator, with the change of 
variables (2.12), leads to a matrix p given by 



/0 1 1 0^ 

1 1 

1 1 
110 

Vo 1 1 0/ 



(5.2) 
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We find that N 5 = 1, and 



A 5 = ^(a? - 2a;a i+ i + 2a;ai+ 2 ) , 



Vij 



Cti-2 - OLi-\ + «i - + a i+2 , 

1 — 2 — 2 <5i,j+i . 



(5.3) 



Plugging the one-mass box integrals (4.36) into equation (3.24), and using the dilogarithm iden- 
tity (4.38), we obtain 



.7 = 1 



7T 

IT 



(5.4) 



In terms of momentum invariants, the unreduced integral is 
rr (-S5i) e (-si2) e 



h = 



(- S 23) 1 + e (-S34) 1 + e (-S45) 1 + e U 



^ + 2Li 2 (l-^) + 2Li 2 (l-^)-?L 

e^ V s 5 i/ V si 2 / 6 



cyclic + 0(e). 
(5.5) 

From this expression we can obtain the value in any region by using the usual ie prescription and 
observing that 1$ is manifestly real in the region where all Sij < 0. 

For the tensor integrals, we do need some information about the 0(e) parts of the scalar 
pentagon. It turns out that leaving the six-dimensional pentagon j^= 6 - 2e [ n equation (5.1) leaves 
us with enough information about these terms that we can use the scalar pentagon as a generating 
function for the tensor integrals to 0(1), without having to evaluate /^ ,=6_2e explicitly. (The 
explicit solution for I® =& involves a rather long combination of Li 3 's, Li 2 's, and logarithms whose 
arguments are complicated solutions of various quadratic equations.) We show how to do so in the 
next section. 

5.2 The Pentagon Integral with One External Mass 

For the pentagon with one external mass, m 5 / 0, we use the same change of variables (2.12), 
(2.13) as in the massless case, except that we also define the rescaled mass m| = — a^aiml, which is 
taken to be a variable independent of the . We find that the normalization factor is A^ 5 = 1 - 
while the rescaled Gram determinant is given by 



- m 



5' 



Ag = A^ m + rng(-2aia 3 + 2a 2 a 3 - 2ai, - 4a 2 a 4 + 2a 3 a 4 - 2a 3 a 5 ) + QjK) , (5.6) 

where Ag m is given in equation (5.3). Using these values in the general expression for the scalar 
pentagon (3.24), and collecting the box integrals with one and two external masses from section 4, 
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we get 



4 m W = rr(n^)l? 

+ 2a 2 



2„, \„2 e l 1 ~ (™5) 



a\ +2e + a^ +2e + a\ +2e + [(«i - m^ 2 )a? £ + (a 5 - m§a4)a§ e ] — 



2\-e 



Li 2 1 



+ 2a 3 



«2 \ 



+ Li 2 1 



«3 

a 2 

Li 2 I I - ] + Li, ( 1 — — 



+ 2a 4 



Li 2 1 



«5 

Q?4 



+ Li 2 1 



«3 
«4 



«3 / 



a 3 



«i - a 2 + (1 - mg)Q 3 - a 4 + a 5 , „ 2 \ 
+ : r-o Li 2 {l-m£) 



+ 2 



cci — m 5 a 2 
1 — m| 



a 2 



Li 2 1 Li 2 1 m 



"2 . 2 



ai 



+ 2 



1 - m 2 
\ — m\ 



Li,[l-^ 

a 5 



7T 



a 3 } + 0(e), 



Li 2 I 1 m 5 

a 5 



(5.7) 



or in terms of more conventional kinematic variables 



4 m [l] = - 



S45S51 



-S34 



-S45 



+ 



(-545) 1+e (-55l) 1+e , (-^5l) 1+e (-Sl 2 ) 



1+e 



+ 



(-S5l) e (-Sl2) e (-S 23 ) e (-Sl2) e (-S 23 ) e (-S34) e (-^("^H"^ 



+ 



S45S51 - mgs 23 



S23S34 1 



mt 



S51 



(-S23) e (-^34) e 



-S45J 



-Sl2j 



(-Sl2) e (-«23) 



+ S12S23 1 If , e( w X 

S45/ (-S34) £ (- S 45) e (-S5l) 



1 _ /^|f23\ 
VS45S51/ 



+ 2s 34 s 45 



S23 



+ 2S45S51 



Li 2 1 



S45 
S34 
S51 



Li 2 1 - — + Li 2 1 - — 



S51 



+ Li 2 1 



S34 
£12 
S45 



+ 2s 5 iSi2 



S23 



S51 



Li 2 I 1 — — ) + Li 2 I 1 — — 



S12 



+ S45S51 



■523 s 34 — "S34S45 + "S45 s 51 ~~ m 5 S 23 ~ *51 s 12 + Sl 2 S23 



Lio 1 



S45S51 - "I-5S23 



ml 



£23_\ 



S45S5I / 



+ 2 
+ 2 



(ggl - m^)s 2 3-S34g45 

S45S51 - S23"T-5 
(s 45 - m|)£5iSi2£23 

S45S51 - s 2 3ml 



Lio 1 



S45 
S23 



Lio 1 



Lio 1 



^]-Li 2 

S23 



mt 



S51 
S45 



7T 



y s 45S 5 i > + 0(e). 



(5.8) 

Observe that it has the expected symmetry under flipping external legs 1 4 and 2^3. The limit 
of the expression (5.8) as 777,5 —> does not yield the massless pentagon integral (5.5), for similar 
reasons as explained at the end of section 4 for box integrals. The single mass pentagon I^ m , which 
is given through O(e ), should not be expected to have a smooth limit onto the massless pentagon 
as 7775 — > 0, because Ig m incorporates the box integral I% mh , and through it the triangle integral 
/| m which does not have a smooth limit. 
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6. Feynman Parameters in the Numerator 

In this section, we explain how to use the scalar pentagon I5, when expressed in terms of box 
integrals and jjp- e ~ 2e v ia equation (5.1), as a generating function for the tensor integrals I^[P(ai)] 
through O(e ). The general discussion applies to the pentagon integral with any number of external 
(or internal) masses; we shall also give explicit formulae for the massless pentagon at the end of the 
section. 

The only complication in applying the differentiation formula (2.19) is the appearance of 
I,P =6_2e and its derivatives at O(e ) when the degree of P(a,j) is two or higher. It is easy to 
eliminate the derivatives of j^— 6 ~ 2e [ n favor of j^— 6 ~ 2e itself and the D = 6 — 2e scalar box 
integrals I^~ 6 ~ 2t ^\ We just let e — > e — 1 in equation (3.17), whence 



ofD=6-2e 

= (-1+&) 



da 



1 



1 („ 7i7j \ fD=6-2e (j) li f D =6-2e 



(6.1) 



Since each term in this equation is nonsingular as e — ► 0, and since we need if 6 2e only to 0(e°), 
we can set e = in jf= e ~ 2e anc l if~ 6 ~ 2e ^\ and use in place of (6.1) the slightly simpler equation 

^P- = -T— (na - ^) it 6 (J) - 2- 4 D=6 • (6.2) 

The D = 6 scalar box integrals can be worked out directly, or they can be determined from 
the D = 4 — 2e box integrals and triangle integrals, using equation (3.25) with n = 4. For the box 
with one external mass, needed for the massless pentagon, the explicit result is 

It" U) = "-r 1 ^ , (6-3) 

A5-7J 

where Lj is defined in equation (4.43). 

Having eliminated its derivatives, we still have to deal with the appearance of i 5 D=6 itself in the 
integrals Is[P(ai)], for m > 2. The way to proceed is suggested by an argument due to Ellis, Giele 
and Yehudai [23] . They work in terms of loop- momentum integrals directly, and use the Brown- 
Feynman or Passarino-Veltman procedure to solve for the tensor pentagon integrals in terms of 
lower-order tensor integrals (pentagons and boxes), all evaluated in D = 4 — 2e. The quantity 
I® =6 does not appear at O(e ) in any momentum-space tensor integral. This fact suggests that 
in our approach, 1 5 D=6 will cancel out of the integral of any Feynman parameter polynomial that 
is the Feynman parametrization of some tensor integral in momentum-space. In appendix III, we 
show explicitly that this is indeed true for integrals with up to three loop-momenta inserted. It is 
straightforward to extend the argument to five loop-momenta, the maximum number encountered in 
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any gauge theory amplitude. (Beyond five insertions of the loop momentum, ultraviolet divergences 
of the integrals complicate matters.) 

While /f =6 disappears from the final answer in any gauge theory calculation, it is still useful 
to know with what coefficient it appears in any particular term. One may use the vanishing of its 
coefficient in the final expression as a check on the complete calculation. Also, it is simple to write 
recursive formulae for the integrals of monomials in the Feynman parameters using this information. 
Let us work out the coefficient of I® =& in I^a^ . . . a» m ]. Define di 1 ...i m by 



I 5 [a h ...a im ] = 



- f D = 6 + 

_1 J 5 > ■ ■ ■ i 



(6.4) 



N 5 A™ 

where '. . .' denotes scalar box integrals (in D = 4 — 2e and in D = 6) and their derivatives. Notice 
from equations (6.2), (3.8), and (3.10) that \f^i^ =& satisfies a simple equation, 



da; 



= --Y. 



i=i 



Vij 



'A 5 7f= 6 



Now write the term (e/N 5 ) A 5 I® 6 in equation (5.1) as 



e 

^5 



and apply the differentiation formula (2.19) to get 



(-l)-A 



m A m_1 / 2 
5 



1/2 



m > 2. 



(6.5) 



(6.6) 



(6.7) 



2(m - 2)! da^ . . .da im 
We have taken the limit e — > in the T-function prefactor in (2.19), since we are working only to 
O(e ). Carrying out the differentiations explicitly for the cases of interest, m = 2,3,4,5, we get 



d - = ± 



Hjk 



j 

u ijkl " 



d 



ijklm 



Vij A 5 - lilj 

(Vij Ik + Vjk li + Vki lj) A 5 - 3 liljlk 

■ {VijVkl + VikVji + mv.jk)Al 

- 3(r/jj7fc7/ + Vikljli + ViHjlk + Vjkiai + Vjililk + mnaj)^5 + 157i7j7fc7i 

■ {vijVknm + perms of ijklm (15 terms)) Ag 

- 5{vij7kin m + perms of zjfcZm (10 terms)) A 5 + 35 7^7^7/0 7z7« 



(6.8) 

In some calculational schemes for gauge theory amplitudes, the D = 6 — 2e pentagon integral 
will itself appear with Feynman parameter polynomials of degree m < 3 inserted. It is then useful 
to know the coefficient df =< ? defined by 



?D=6 



L>=6 



1 %\...%r, 

Oj, . . . Oj = — ^ — 

n lmJ Am 

^5 



(6.9) 
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where again '. . .' denotes scalar box integrals (in D = 6) and their derivatives. Writing />P~ 6 = 
A^" 1 ^ 2 x (A^ 2 I® =6 ) and repeating the above steps we find 



d 



D=6 



(-l) m A 



m+l/2 
5 



d m A~ 



-1/2 



" i "" m ml dcti 

The explicit values for the cases of interest are 

di = 7i , 



m > 2. 



(6.10) 



d 



ijk 



J _ _1 

^7. 7 2 



Vij A 5 - 37i7j-J , (6.H) 

7fc + ??jfe 7i + »7fci 7j) A 5 - 5 7i7j7fc • 

Now we shall give explicit formulas for the massless pentagon integrals with up to two Feynman 
parameters inserted, along with a simple recursion relation for generating the remainder of the 
integrals. 

For a single parameter insertion, equation (2.19) gives 

1 dh 



1 + 2e dcti 

Thus we may differentiate the O(e ) expression (5.4) for j 5 , using also equation (4.37), to get 

CXi-l 



(6.12) 



h [at] = r r 



% + 2Li 2 fl-^±i 



+ 2 Li 2 1 - 



OLi 



7T 
~6 



2 1 



+ 0(e). 



In the case of two Feynman parameter insertions, we have 

f r l 1 92/5 



(6.13) 



(6.14) 



2e(l + 2e) 9a, da 3 ' 

which must now be applied to the expression (5.1) for I§. The 0(1) terms in 1$ [a^aj] receive 
contributions both from I§ =& and from the 0(e) terms in the box integrals Since 7, is linear in 
the cti, only derivatives of the 1^ appear on the right-hand side of (6.14). The derivatives dljf^ /dcti 
at 0(e) are nothing but the single insertions /^[aj] at 0(1), thanks to the box differentiation 
formula 



if [ 



1 31 



(?) 



(6.15) 



4 L ~ lJ 2e d ai 

These integrals are tabulated in equation (4.42). Carrying out the differentiations in (6.14), we find 

f, 1 / kj f a i+i a t-i . a\_ 2 {al - a\ +l ) a? +2 (a| - af^) 
hyaidj] = rr<—g-[— 1 T 1 T 



5 i+ ij a\_ 2 + 8 j+1 ,i a £ j_ 2 ct\ - a) 



+ E 



k=i 



1]ik7j + Vjk^fi ~ VikVjklk 



liljlk 



d 



+ *L j£=6 J, + (e). 



As - it A 



(6.16) 
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For a generalization of this formula to arbitrary pentagon kinematics, and also to hexagon (n = 6) 
integrals, see equation (VII. 8) in appendix VII. 

For more than two parameters inserted, we can proceed recursively. Define some new quantities 



jphys 



cii 1 . . . a>i m ] to be the "non-/,p 6 " terms in I^a^ . . . Oi m ], i.e. 



I5 [ a ii ■ • • a i 



jphys j 



+ 



d 



A' 



m— 1 5 



(6.17) 



Then the differentiation formula (2.19) along with (6.2) generates the following recursion relation 



fphys r 



-1 



m - 2 - 2e 



dai 



+ 4- 



i=i 



~ -7? 



(6.18) 



In applying this formula, it is convenient to have a differentiation formula for the Lj, in terms of 
logarithms: 

a j+1 ln(a j+1 /a j+2 ) 



dai 



a, 



(k. 



3+2} 



3-2) 



a j+1 - a j+2 
aj-i ln(aj_i/a J _ 2 ) 



+ \n(ctj-i/aj-2) 



+ Ha j+1 /a j+2 ) 



(6.19) 



Otj-l — CHj-2 

This completes our prescription for evaluating massless pentagon integrals with Feynman pa- 
rameters inserted, in terms of dilogarithms and logarithms. The same basic procedure also works 
when external and/or internal masses are present, provided that the relevant box and triangle in- 
tegrals are known through O(e ). (The triangles appear through equation (6.15) in combination 
with (3.16).) If all internal lines are massless, then all the requisite boxes and triangles can be 
found in section 4, except for the three-mass triangle. This triangle may be computed in D = 4; 
see for example refs. [6,17]. In appendix V it is computed in D = 4 — 2e for arbitrary e, as a further 
illustration of the partial differential equation approach to scalar integrals. 
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Appendix I. Collection of Massless Pentagon and Scalar Box Results 

In this appendix we collect those results that are useful in an explicit calculation. The massless 
pentagon integral of interest is 



h[P m (M)] = T(3 + e) / d 5 

Jo 



5(l-Zi^)Pm({ai}) 



-si 2 aia 3 - 52302^4 - 5340305 - s 45 a 4 ai - s 51 a 5 a 2 - ie 



3+e 



(LI) 
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where P m ({ai}) is a polynomial in the aj of degree m. For use in differentiation formulae we define 
the reduced integrals 

i n [P({ ai })] = m«,) I n [P({ ai / ai })] . (1.2) 
where Sj^+i = — 1/(0^0^+1) mod 5. The basic differentiation formula for the pentagon is given by 



P m ({ai}) 



r(2-m + 2e) 



T(2 + 2e) 

Through 0(1) the scalar pentagon is given by 
5 



dotj 



h [1] • 



(1.3) 



.7 = 1 L J 



+ 0(6) 



(1.4) 



where rp = T(l + e)T 2 (l — e)/r(l — 2e). The case of one Feynman parameter in the numerator 
may be obtained by directly applying the differentiation formula (1.3). 

Beyond one Feynman parameter it is best to use the explicit value of the two Feynman pa- 
rameter integral as a generating function for integrals with three or more Feynman parameters in 
the numerator. The two parameter integral is given by 



h[a,iaj 



Oti 



ati - a i+1 



cti - ati-i 



S i+l,j a i-2 + Ol)-2 a i - a j 



on — a-: 



+ E 



k=i 



Vikjj + rjjkli ~ VikVjklk 



liljlk 



A, 



Li 



5 J A 5 - 7fc 



A 5 



+ ^ I 5 D = 6 



+ 0(e), 
(1.5) 



where 



A 5 = ^2(a 2 j - 2 aj a j+1 + 2a j a j+2 ) 
15A 5 



7i = o 



Vij = 



dij 



2 dat 
chn 1 d 2 A 5 
daj 2 daidaj 

Vij A 5 - lilj 



- cti-i + ai - a i+ i + a i+2 , 



(1.6) 



-1, * = J±1, 
+1, otherwise, 



and 



= Li 2 1 - 



«i+2 



+ Li 2 1 - 



ai-2 



+ ln 



«i+2 



In 



CCi-2 



7T 

~6 



(1.7) 



In calculations there is no need to know the explicit value of the six-dimensional pentagon 
Lf =6 since it cancels from all quantities arising from loop momentum integrals. However, when 
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applying the differentiation formula (1.3) the terms containing I® 6 cannot be dropped since they 
generate D = 6 box integrals via the equation 



da. 



ST 1 fr> 7i7 ^ ^ f D = 6 7i f D = 6 



(1.8) 



Useful formulas when applying the differentiation formula (1.3) are 



9L± 



+ 



1 

aii 
1 



{Sij-i - Sij-2) 



a j+1 ln(a j+1 /a j+2 ) 
aij+i - a j+ 2 

atj-i ln(q J _i/a J _ 2 ) 
oij-i - a.j-2 



+ \n(a j - 1 /a j - 2 ) 
+ Ha j+1 /a j+2 ) 



(1.9) 



and 



?D=6 (j) 

-1 A 



4 



A5-7? 



(1.10) 



where Lj is defined in equation (1.7). 

We collect here the dimensionally-regulated scalar box integrals with massless internal lines, 
but 0, 1, 2 or 3 nonzero external masses, which appear in the process of evaluating (n > 5)- 
point integrals, and in subdiagrams in QCD loop calculations. The integrals are defined through 
equations (2.5) and (2.6). 
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Appendix II. Connection with the Work of van Neerven and Vermaseren 

Melrose [5] and Van Neerven and Vermaseren [7] were able to represent the general scalar 
pentagon integral in D = 4 as a sum of five D = 4 box integrals. On the other hand, equation (3.23) 
expresses the pentagon integral in D = 4 — 2e as a linear combination of five box integrals (also 
in D = 4 — 2e), plus the pentagon in D = 6 — 2e dimensions, so it can be thought of as the 
dimensionally-regulated version of the equations in refs. [5,7]. Indeed, the D = 4 equation in 
ref. [7] was our motivation to find an algebraic D = 4 — 2e equation. (Similar relations have 
recently been found using momentum-space, rather than Feynman parameter, techniques by Ellis, 
Giele and Yehudai [9].) We would like to verify that the D = 4 — 2e and D = 4 equations 
are consistent with each other, or in other words that the (D-independent) coefficients of the 
box integrals in equation (3.24) are equal to the corresponding coefficients in refs. [5,7] (up to 
normalization conventions for the integrals). To do this, it is simplest to rewrite equation (3.24) in 
terms of unreduced integrals as 



h[l] = 5 /<*>[!] + 0(e) (11.1) 



i=l 

where 



Ci 



The second form of the Cj, in terms of more conventional kinematic variables (the matrix S is 
defined in (2.4)), is the form in which the Cj were obtained in ref. [10]. In this form the are 
manifestly the same as those found by Melrose. 

The coefficients found by van Neerven and Vermaseren involve the D = 4 Levi-Civita tensor, 
and are not manifestly equal to (II. 2). Expressed in our notation, they are given by 

4A 5 - 2 J2t=i v i' w 2vi-w . 

Cl = vfl - 4A 5 M? ' C * +1 = ~ w 2 — 4A 5 Mf ' * = 1 ' 2 ' 3 ' 4 ' ( IL3 ) 
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Here the "axial vectors" Vi are the D = 4 duals of the vectors pi appearing in the momentum-space 
version of the pentagon integral 

yt 1 = £ W2P3P4 yV = £ PlHP3P4 yV = £ PlP2HPi yL 1 = £ PlP2P3^ 

1 ^2 '3 ^4 ^ 

- E fc ;> p" = o. (IL4) 

where £ W2P3PA is short for ^n^^mP^P^P^ ■> etc The Gram determinant of the vectors pi is 
A 5 = e PiP2P 3 P4 £pip2p3p4) and i s defined by 

4 

^ = J><, = p 2 + M 2 -M 2 +1 , i = 1,2,3,4, (II.5) 

i=l 

where Mj are the masses on the internal lines. The definition of the Cj in (II. 3) may seem to be 
tied to D = 4, because of the presence of the axial vectors. However, the inner products Vi ■ Vj can 
be eliminated in favor of the inverse of the matrix tij = 2pi • pj , according to 

Vi ■ vj = 2A 5 (t-% , i,j = 1,2,3,4. (II.6) 

Thus the q can be written in a D-independent form: 

"2 + 2J2l i=1 (t-% rj 



Cl 



4 (t-i)-r- (IL7) 

Ci+i = —7 — , i = 1,2,3,4. 

EUiM*- 1 )*,!-, - 2Mf 

To show that the Cj in (II. 7) agree with those in (II. 2), it suffices to show that they obey 



Y^SijCj = 1, i = l,...,5, (II.8) 

since 5 is generically invertible for n = 5. When internal masses are also present, S is given by 

S l3 = \{Mf + M 2 - (pi_! - Pi -i) 2 ) = i(2M x 2 - ri _i - r,_! + ^-ij-i) , (II.9) 

and Vj-i = ti-ij-i = for j = 1, so that 

5 1 , 5 . 1 4 

E% c ^ = 2 (2M ' -^(E^) + 2^ ( "^ +tij)c '' +i) • (IL10) 

j=i S=i y j=i 

Plugging in the values of q from equation (II. 7), we get 

(~2Mf + r^i) + (£ 3 4 >fc=1 r, {t~ l ) ]k r k - n . x ) 
H ^ ' " EUirikCt-^wr, - 2Mf 

as required. 

In the same fashion, an equation obtained by van Neerven and Vermaseren, relating hexagon 
integrals to pentagon integrals, can be shown to be equivalent to equation (3.23) for n = 6 (and 
D = 4 external kinematics). 
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Appendix III. Proof That if" 6 Drops Out 

An explanation of why an explicit computation of I^ =6 is not needed for the evaluation of 
pentagon integrals near D = 4 comes from the momentum-space representation of tensor integrals; 
when performing a Passarino-Veltman decomposition IP =6 never appears [23] and therefore it 
can be expected to cancel from amplitudes evaluated using the Feynman parameter techniques 
discussed in this paper. In this appendix, we will demonstrate that lP =6 cancels when summing 
over contributions which reconstruct the loop momentum integrals appearing in dimensionally- 
regulated four-dimensional field theory amplitudes. Thus, there is no need to explicitly evaluate 
lP =6 . (In this appendix we treat lP =6 as equivalent to I^~ 6 ~ 2e ; since lP =e is completely finite 
the difference between the two is of 0(e).) The argument holds for general kinematics (arbitrary 
external or internal masses), though here we suppress internal masses. 

Define the general pentagon integral by 



J 5 [P(p")] 3(4^) 2 -4! J -±^L J d 5 a 



(aip 2 + a 2 (p-pi) 2 + a 3 (p-p 2 ) 2 + a 4 (p - p 3 ) 2 + a 5 (p-p 4 ) 2 )' 

(III.l) 



where pi = 5^}=i kj an< ^ P{p^) is some polynomial in the loop momentum p^. The normalization 
factor in front ensures that the integral, when Feynman-parameterized, is normalized in the same 
way as the integrals I 5 [P(a,i)] defined in section 2. 

In order to relate the integral (III.l) to Feynman-parametrized integrals of the form (2.5), we 
complete the square and integrate out the loop momentum in the usual fashion. To complete the 
square in the denominator, we shift the loop-momentum variables to 

4 

P = <? + E a-i+iPi ■ (III.2) 
i=i 

Integrating out the loop momentum, for up to three powers of loop momentum in the numerator, 
then gives 

4 

h\p»]=Y,h[a i+1 } P ?, 

i=i 

i 4 

i i i 

+ ^PiP v 3 P P k h[ai+ia J+1 a k+1 \ . 

ijk 

(III.3) 
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Here we will explicitly consider only up to three loop momenta; the other cases follow similarly. 

For the case with one loop momentum inserted, since the explicit value for I 5 [a i+ i] given in 
equation (6.13) does not contain I® =6 , there is nothing to check. Beyond this, we have from 
section 6 that the coefficient of />f =6 in the explicit value for I^[a\ ■ ■ ■ afc] is given by 

(-l) m A5 /2 d m Al /2 



«^ 5- a „ & a „ (m > 2), (IIL4) 



-.i. a -. ra 2iV 5 (m-2)! 11 12 lm d ail da i2 ■■•da i 
so that for m = 2, 3, 

c i.?fc = ? „ Wj 7fc + %fc 7i + »7fei 7j)A 5 - 3 7i7j7fc 



''' 2iV 5 A| 



(III.5) 



where 7 j and rjij are defined in equations (3.8) and (3.10). 

The identity that we will use to show that I^ =6 cancels is 



2N 5 A 5 



raj A 5 - 7i7 ,) = (t-V^-i . *,J =2,3,4,5, (III.6) 



where t^- = 2pi • pj. To verify the identity, we multiply it on the right by ij— i, which can be 
written [10] in terms of the matrix p = N n ri~ l using 

Pi-i-p,--! = ^_^_^ + ^L, z,j = 2,3,...,n. (III.7) 

ttjCKj a^ai ct\a.j a\ 

Thus we have (using the equations (3.9)-(3.13) that relate p, rj, 7 j and on) 

5 5 

E \ - aiaj / - _ \ f _Pjk_ Pji Pik . PiA 
.J* 1 *- 1 *- 1 2- ,v,A, V''" 7 * 7 vUi«* a?J (HI.8) 

= (5 ifc , i, fe = 2,3,4,5. 
Equation (III. 6) implies that 

£ P^+U+i = E = ^j? > ( IIL9 ) 

since the four vectors span D = 4 Minkowski space. 

Using this identity and keeping only the ijf =6 content we then have 

4 



= -^5 D=6 [l] ^4-2e] + E P?P-C i+ l, j+ lI°= 6 [l] + boxes 
1 

2< 

boxes + 0(e), 



= --<5[_ 2e] / 5 D=6 [l]+ boxes + 0(e) 



(111.10) 
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so that I® =6 drops out as claimed. To arrive at the last line, we used the finiteness of I^ =6 and 
that <5p 2e ] can yield only 0{e) contributions. This shows that there are no 'left-over' pieces of I® =& 
remaining when all the pieces are combined to form an amplitude derived from a loop momentum 
integral with up to two powers of momenta in the numerator. 

The three Feynman parameter case is similar. Again applying the identity (III. 6) we have 

/5(PW) = ~\ (5-2e] J>f/ 5 D=6 [a l+1 ] + cyclic) + £pfp>£c^/f= 6 [l] + boxes 

i ijk 

= 4fe E^^x^^w + c H + boxes + °^ (IIL11) 

i 5 

= boxes + 0(e), 

where we used 

Cijk = <kj + cyclic (111.12) 

and 

/f= 6 [ ai ] = ^2i/f= 6 [l]+ boxes, (111.13) 

from equation (6.11). 

It is straightforward to continue in this way, demonstrating that I® =6 drops out from the loop 
momentum integrals encountered in relativistic field theories. For gauge theories, up to five factors 
of the loop momentum in the numerator can appear. 



Appendix IV. Constants of Integration for Box Integrals 

In this appendix we evaluate the constant of integration for the box with two adjacent massive 
legs, or with three massive legs, by performing the integral at the point where all the on are equal. 
The constant of integration for the adjacent two-mass box is a special case of that for the three-mass 
box, with A = 0. We have 

Io ee 7f>, = 1,A) = r(2 + e ) f d* Ul . (iv.l) 

J [{Ui + \U 2 )U 3 + U4{1 - U4)\ 

We let 

ui = z(l-y), u 2 = (1 -z)(l -y), u 3 = y(l-x), u 4 = xy. (IV. 2) 

The z integral is elementary and leads to 

f ° = ~~iZT\-£ dx f ^T=9 { [(l-x)(l-y)+x(l-xy)] - 1_e - [X(l-x)(l-y)+x(l-xy)] "^j . 

(IV.3) 
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The y integral can be done in terms of hyper geometric functions, 

^f^/^u+^i-.i-^) (iy4) 

For A / 0, the integrand has no singularities as e — ► 0, so we may expand it in e; the 
hypergeometric functions have the following expansion for small e, 

2 Fi(l + e,-e;l-e;«) = 1 + e ln(l - v) + e 2 [-2 Li 2 (i>) - \ ln 2 (l - v)] + 0(e 3 ) (IV.5) 
(we only need the first two terms here), which leads to 

I = r r dx 
Jo 



' -±-lnx-ln(l-x) + 21n(A+(l-A)x) 21n(A + (1 - A)x) 



A + (1 - X)x 



+ 



(l-A)(l-x) 



rr ^ + / 1 ^(-ln(n-A)-ln(l- U )+21n U + 21n(l-A)) + 2 £ A ^ln(l-t;)| 



1 - A I e 

rr 



A U 



1-A 



^ - ±ln 2 A 
e 2 



(IV.6) 



In the case A = 0, we add and subtract terms in (IV. 4) to obtain 
T(l + e) f 1 dx 



Jo (A = 0) 



1 — x 



2 Ji(l + e, — e; 1 — e; 1 — x + x ) — x e 2 Ji (1 + e, -e; 1 - e; x) 



F(1 + e) / x" 1 " 6 ( 2^(1 + e, -e; 1 - e; a) - 1) - ^ + ^ / 
e Jo e Jo 



l-e 



(IV.7) 

In the first and second integrals, the integrand is again nonsingular everywhere, and we can expand 
in e; the third is elementary: 



7 (A = 0) = 2T(l + e) / dx 



lnx 

1 — x 



r(l + e) / dx 



+ 



r(l + e)^ vr 2 2 



(IV.8) 



*2 



so that cq(0) = 0. 



Appendix V. The Triangle with Three External Masses 

The differential equations approach also provides an easy way to derive a compact expression 
for the three-mass triangle integral to all orders in e. (The integral is in fact finite, so only the 
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leading order is needed in practical calculations; but in order to examine explicitly the limit in 
which one of the external masses vanishes, it is convenient to have the forms derived here, or ones 
equivalent to them [24].) 

The three-mass triangle with massless internal lines satisfies the following system of equations 
(using N 3 = 1): 
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Also, the two-point integrals are very simple, 
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dot; 



= a; 



-1/2 



7i7j 
A3 



a i+ iai_i = i(7i + 7i_i)(7i + 7 i+ i) = ±A 3 (1 + <$?). 
Therefore we may solve the differential equations (V.l) by 



(V.l) 



(V.2) 



(V.3) 



(V.4) 



(V.5) 
(V.6) 



h = A" 1/2+e [F(8 1 ) + F(5 2 ) + F(5 3 ) + C] , 



(V.7) 



where F(S) satisfies 



F>(8) = -I(l- 2e )Aj 



e(l - 2e) 



(iA 3 (l + 0) 



e-l 



2 l-2e 



(l + <5 2 ) , (V.8) 



and C is a constant of integration. 
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We need the integral 



e-1 



/ dzil + z 2 ) 6 - 1 = / dz (l + iz) e - 1 (l-iz) e - 1 
Jo Jo 

dw w e - 1 2 £ - 1 (l-w/2) 

(l + iS) e 2 F 1 (l-e,e;l + e 
1 + iS 



2 e ~ 1 i r 



e 



1 + iS 



4^-1 i 



e a 



l-i<5 

1 + z<5 
1 -iS 

n-i8\ 



F 1 2e, e; 1 + e; - 



2 Fi ( 2e, e; 1 + e 



J J 1 *-) 



2 

l + t<5 

i — i<y 

1 -iS 
l-iS\ 
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(V.9) 



where we have symmetrized the result in the last line. Alternatively, we may do the integral as 



[ dz (1 
Jo 
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2\e-l 
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OO 
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(V.10) 



(§)■■■ (m-£)(m+£) m! 



777=0 

OO 

m=0 

= <5 2^(1-6, i;f;-5 2 ). 

The two expressions for the integral can be related using a variety of hyper geometric identities. 
Thus we have 



h{«i) = -\^k- l/2 ^[f{5 1 ) + f{5 2 ) + f{5,)+c] , 



(V.ll) 



where 



f(5) = etf-'S 2 Fi(l-a ; |;-<i 2 ; 



1 



2 F X 2e,e;l + e; 



1-H<5\ /l — i<5\ 



l-i5\ 



(V.12) 

To fix the constant c, it is easiest to consider the integral at the following, somewhat asym- 
metric, kinematic point: 



-11 -11 -1 

Sl2 = = -„ , S 23 = = -- , S31 = 

0301 2 ol\ol 2 2 a 2 a.3 



(V.13) 
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or 

Qi = 2, a 2 = «3 = 1; A 3 = 4; <5i = 0, <5 2 = <5 3 
At this kinematic point, we make the change of variables 

ai = 1 - y, a 2 = xy, a 3 = (1 - 



with Jacobian equal to y, and obtain 

4(2,1,1) = ^i±^ f ,1 



= 2T(l + e) / dx / y" £ [l - (1 - 2x(l - " 1_£ 
= 2 e F( | 1 + ^ / dx 2 i ? i(l + e,l-e;2-e;l-2x(l-x)) . 



; r(i + e) 

'0 

Next we use the change-of- variables, 

1 - 

X = 2 

and a hyper geometric identity, to get 
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We then use the following hyper geometric identities, 
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r(l- e )r(i; 



to get 



4(2,1,1) = - 2 l£ 2j F 1 (l- e l;f;-l) + 



r(l-e) 
4 £ vrr(l + e) 



e - • 2e 

On the other hand, plugging the values of A 3 and Si from (V.14) into equation (V.ll), 
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Comparing equations (V.20) and (V.21), we find 

2vrer(l + e) 



= -2vre 



r(l - 2e) 



r r T 2 (l - e) ' 

Despite its appearance, equation (V.ll) does have a finite limit as e — ► 0, 

3 



E 

3 j = l l 



Li 2 



1 + i<5j 
1 - iSj 



- Li 2 - 



1 + iSi 



+ O(e), 



(V.22) 



(V.23) 



which is the form given in ref. [17]. 



Appendix VI. Higher-Point Scalar Integrals 

In this appendix, we discuss formulae allowing the evaluation of higher-point scalar integrals 
(n > 5), in part to correct some statements we made in a previous paper [10]. The corrected results 
will be similar to results obtained previously by Melrose, and by van Neerven and Vermaseren [5,7]. 
The main difference is that the present results allow for external kinematics in the full 4 — 2e 
dimensions, which is useful for obtaining tensor integrals by the differentiation method discussed 
in sections 2, 6, and appendix VII. 

We begin by recalling equation (3.23), which we rewrite here in a slightly different form, 

n i 

f.-Ejfce + f-IH-*)^**. (VLl) 

For n > 6, in order to use equation (VI. 1) to evaluate scalar integrals, it is desirable to take the 
external momenta k±, k?, . . . , k n to be restricted to D = 4. The loop momenta have to remain in 
D = 4 — 2e in order to regulate infrared divergences. In the 't Hooft-Veltman variant of dimensional 
regularization, the external momenta appearing in the one-loop integral in a next-to-leading-order 
calculation are indeed taken to be four-dimensional. In the conventional dimensional regularization 
scheme, the external momenta are taken to be 4 — 2e-dimensional, but this will generally lead to 
only 0(e) corrections, since the integrals I® =6 ~ 2e are finite as e — ► for n > 4. 

In reference [10], we argued that the term containing I^ =6 ~ 2e in equation (VI. 1) could be 
dropped for n > 6, when the external momenta are restricted to D = 4. The argument was based 
on the fact that for n > 6 the Gram determinant A n appearing in equation (VI. 1) vanishes for 
D = 4 kinematics, due to the linear dependence of the (n— 1) vectors k\, &2, • • • , k n -\ [18,25]. If the 
I^ =6 ~ 2e term could be dropped, then equation (VI. 1) would reduce to a simple recursion relation 
expressing the scalar integrals /„ as a linear combination of the n (n— l)-point integrals I n —\. For 
n = 6, the argument does indeed hold, and the scalar hexagon integral is given by 

6 

J 6 = ^27^5° (D = 4 kinematics). (VI.2) 

i=i ^ 
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Unfortunately, for n > 7 the situation is more complicated. It is true that for n > 6 the Gram 
determinant A n vanishes for D = 4 external kinematics. However, for n > 7, the factor of N n in 
the denominator also vanishes. Indeed, N n is given by N n = 2™" 1 detp = (Uti "0 det5 > 
and the dimension of the null space of the n x n matrix Sij is n — 6 for D = 4 kinematics [5]. 
Therefore, for n > 7, the coefficients appearing in equation (VI. 1) are not well-defined for D = 4 
external kinematics (which is where we would like to use the equation). 

Notice that both numerator and denominator of the coefficient ratios ^i/2N n and A n /2N n 
vanish for D = 4 kinematics: The matrices that give rise to A n and to 7$ = ^(dA n /dai) have null 
spaces of dimension n — 5 and n — 6 respectively. Based on the dimensions of the corresponding 
null spaces, we can argue that A n vanishes "faster" than N n , and 7, vanishes "equally fast", as 
D = 4 kinematics are approached. Thus we might expect that a modification of equation (VI. 1) 
should exist, which is well-defined for D = 4 kinematics, and for which the coefficient of I® =e ~ 2e 
vanishes in this limit. In fact, van Neerven and Vermaseren [7] have shown how to obtain such an 
equation, which expresses an ra-point scalar integral in terms of six (n — l)-point integrals. (Their 
derivation was carried out for D = 4 loop momenta; however it is easy to see that it is equally valid 
for D = 4 — 2e loop momenta as well, as long as the external momenta are restricted to D = 4.) 

Here we will obtain an equation similar to (VI. 1), but where the coefficients have N^°2i in the 
denominator instead of N n . Since is nonzero for generic D = 4 kinematics, this equation will be 
well-defined for the heptagon integral (n = 7) in D = 4. It reduces to the above-mentioned equation 
of van Neerven and Vermaseren in D = 4, but it is also well-defined away from D = 4, which makes 
it a useful starting point if one wishes to apply the differentiation approach of sections 2 and 6 to 
compute tensor integrals. The reason why restricting external kinematics to D = 4 complicates the 
differentiation approach is that the variables are then subject to various Gram-determinental 
constraints [18,25], which would have to be respected in performing the differentiations. After 
carrying out the differentiations it is permissible, and usually desirable, to restrict the external 
kinematics to D = 4, in order to take advantages of certain simplifications. An example of this 
procedure, for the one-parameter heptagon integrals, is provided in the next appendix. 

To derive the new scalar equation, we first need some general relations between the quantities 
A n , 7j and iV n , which are associated with the integral I n , and the corresponding quantities A^\, 
7^ fc ' ) and N^ l k } 1 associated with the (n — l)-point "daughter" integral As in section 3, we 

choose the variables for the daughter and parent integrals to be the same. We also take the 
kinematics to be general for now, i.e. not restricted to D = 4, so that all quantities are well-defined. 
The necessary relations follow from the observation: If A is a symmetric n x n matrix, and 
is the (n — 1) x (n — 1) matrix formed by crossing out the k th row and k th column of A, then the 
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inverse of B^ can be computed as 



S"5) = AJ-^*L, i + Kj + k. (VI.3) 

,l i A kk 



The proof is simply to multiply equation (VI.3) on the left by {B^) a = A u , and simplify. Note 
also that (j^fj^j vanishes for either i = k or j = k. 

Starting with the expression (3.8) for /S^-i m terms of on and the matrix rf k \ and using 
equations (3.9) and (3.10), we can rewrite as 



a r] K 'a 



T -1 (fc) _ ly n-l , T J Jk) 



7 77 77^77 7 = — —7 p;p;i!j 

n 



0" 



\ Pkk I 



T 7fc 

7 P7 - -=Y 



Using the definitions (3.9) iV n = 2 n " 1 detp, = 2 n ~ 2 det p^, and the fact that det pW = 

Pkk ^et P * s * ne cofactor of the kk element of p, we have 

Vkk -1 2N n-l nn 

w = pkk = -^r~ • (vi.5) 

Using equations (VI. 4), (VI.5) and the relation 7 T p7 = N n A n which follows from equation (3.11), 

I (k) 

we obtain expressions for A n _ 1 and its derivatives with respect to c^: 

Uk) _ Vkk &n ~ 7fc 

n ~ X "" 2N n 

(k) Vkk7i ~ Viklk /yj 

ll 2N n V ' ^ 

(it) VkkVij - VikVkj 
Vlj ~ 2N n 

One can iterate this procedure to get expressions for A^l^, etc., if necessary. 

Now we proceed to derive the new scalar equation which is of use for n = 7. To do this, we 
consider equation (VI. 1), and also the one-parameter equation 

n 

L[a k ] = E|;^ + ("- 5 + 2e )?| ^ 6 - 2e • (vi-7) 

1=1 n n 
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Multiply equation (VI. 7) by j k /i] kk and subtract it from equation (VI. 1), to get 



E 

i=l L 



2N n Vkk ■ 2N n 



Ik 



Vkk 



A, 



it 



2N n 7] kk ■ 2N n 



?D=6-2e 



(VI.8) 



The coefficients in brackets in equation (VI.8) can now be rewritten in terms of (n — l)-point 
quantities using equations (VI. 6). We get 



n (fe) 
i=l /iV n-l 



Ik 



2N 



(fc) 



■in[ofc] + (n-5 + 2e) 



A (fe) 

"n-l jD=6-2e 



2N 



(k) 



(VI.9) 



re— 1 re— 1 

(fc) 

Any value of A; = 1,2, ... ,n may be used in this formula. Note that 7^ =0, so there are only 
n — 1 terms in equation (VI.9). 

For n = 7 and (generic) D = 4 kinematics, we have Nq^ / 0, while = and j k = 0. So 
equation (VI.9) reduces to 



E 



(fc) "6 



(D = 4 kinematics) , 



(VI.10) 



which contains only six hexagons due to the vanishing of . Indeed, the formula can be shown 
to be equivalent to the scalar integral formula of Melrose, and van Neerven and Vermaseren [5,7]. 
For n > 7, equation (VI.9) is still ill-defined. Presumably one could go on to construct equations 
in terms of ^ k ' l \ A^lj, etc. that will be well-defined for n = 8, and so on. This would be useful 
for evaluating the corresponding tensor integrals via differentiation. 



Appendix VII. Higher-Point Tensor Integrals 

In this appendix, we derive formulae allowing the evaluation of tensor integrals for the pentagon 
(n = 5) and hexagon (n = 6) integrals, for arbitrary internal and external masses. We also briefly 
discuss tensor heptagon (n = 7) integrals. The situation regarding tensor integrals is similar to 
the case of the pentagon discussed in section 6 and appendix III. In order to effectively use the 
differentiation approach, one must show two things: First, that the 1/e pole encountered in the 
basic formula (2.20), at the level of n — 3 Feynman parameter insertions in the n-point integral, 
does not present any problems; and second, that the "hard" six-dimensional integrals I® =6 ~ 2e (for 
n > 5) always drop out of any "physical" tensor integral, i.e. any integral which is the Feynman- 
parametrization of some loop-momentum integral. We discuss these issues here to some extent for 
n = 6, 7; presumably both points can be shown to hold for arbitrary n. 

For n = 5 and n = 6, the insertion of a single Feynman parameter can be treated using 
equation (VI. 7). For n = 5, the term containing the D = 6 — 2e scalar integral l,f =6_2e is 
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0(e) and can be ignored. For n = 6, this term is 0(1), and so we would like to show that 
for "physical" one-parameter integrals (linear combinations of the parameters corresponding to 
Feynman-parametrization of some loop-momentum integral), and for D = 4 external kinematics, 
the integral Iq~ 6 ~ 2€ drops out. Feynman parametrization of the loop-momentum integral I n [p^} 
leads to a linear combination of one-parameter integrals, similar to the first equation in (III. 3), 



(VII.l) 



i=2 



i=2 



So we can show that Iq 6 2e drops out by showing that 



i=2 



(VII.2) 



for D = 4 external kinematics. To show that equation (VII.2) holds, it suffices to contract the 
equation with a set of vectors Pj_ 1 that span D = 4 (we can pick any four of j = 2, . . . , 6 for 
nonexceptional momentum configurations). We then use equation (III. 7) to write Pi-\ ■ Pj-i in 
terms of the matrix p = NqT] -1 , and equation (3.11) to simplify the sum: 



£ 

i=2 



i=l 
6 



N 6 



Pij 



Pn 



Ji_( Pji 



a 



N 6 oti 



3 

Pll 

a 4 



Pa 



OLiO,\ OL\OLj 
Nq «i 



plj , Pll 

~T~ 2 



Pll + Pll 



a,- 



(VII.3) 



+ 



Pu 
ai 



But A 6 = while iV 6 / for D = 4 external kinematics, so 7 6 



D=6-2e 



does drop out as desired. 



We turn next to the insertion of two Feynman parameters. The first part of the derivation 
parallels the derivation of the one-parameter equations (3.4) and (3.16) in section 3. We consider 
the integrals 



Jn;i [^fe] 



rl pl—a n -i p 

T(n-3 + e) / da n -i I da n _ 2 ■ ■ ■ 

Jo Jo Jo 



1 — ai — ai - 



-O.n-1 



at 



n— 3+e 



(VII.4) 



a n = l — a 1 —a 2 - 



-a> n -i 



evaluated two different ways, to obtain the set of equations 



E 

3 = 1 



Pji 



Pnj 



I n [a 3 a k ] 



(n) 



a,,, 



+ 



ai 



°nk 



rD=6-2e 



[1] . (VII.5) 
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Solving for I n [aiCij], we get 



i n [a iaj ] = — |- E [Vjt - ^ ) ft-lfai] + Y-inWi] + 
71 e=i ^ ^ n ' ^ n 



2N„ 



Vij 



7i7j \ fD=6-2e 



A, 



In • (VIL6) 



We can rewrite the right-hand-side of equation (VII. 6) in terms of scalar integrals only, with the 
help of equation (3.25): 



1 n ( 



+ 



+ 



1 



ljl£ 



1 



P =i 



n 

E^n-i + (n-5 + 2 e ) 7l / 5 



+ (n-6 + 2e)7 f)et- 2£ 



D=6-2e 



.t=l 



2N„ 



A. 



(VII.7) 



In this equation, 1^-2 is the (n — 2)-point scalar integral obtained from by eliminating 
the p-th propagator. We keep the original kinematic o^-variables defined for I n ; I n _< 1 will be 
independent of ag and a p . The other quantities — ^n-i an< ^ ^ s derivatives — refer to 

the normalization, rescaled Gram determinant, and so on, associated with I n _i- We can eliminate 
from equation (VII.7) in favor of 1^-2 an d ^n-x^'^i an d use equations (VI. 6) to simplify 
things. We get finally 



InldiClj] 



77 i jA re + (w-6 + 2e)7 i 7j ^ D=6 _ 2e 



2V n A n 
n — 6 + 2e 

+ — E 



+ 



4iV2 



1 n 

n e, P =i 



me 



'n-2 



A r 



fD=6-2e (i) 
1 n-l 



(VII.8) 



This formula merits several comments: 



1 ) The expression f^-2 nas no meaning for I = p; however, I / p is enforced automatically by 
the prefactor. 

2 ) For n = 5, and all-massless kinematics, this equation reduces to equation (6.16); notice that 
r]u = 1 in this case, and that we wrote out the I n -2 terms — in this case triangles — more 
explicitly there. 

3 ) For n = 5 and general kinematics, we now have JsfojOj] to 0(1), which means that we have 
surmounted the "1/e barrier" for the pentagon. That is, insertions of more than two Feynman 
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parameters can be obtained using just the differentiation formula (2.20), and equation (VII. 8) 
evaluated to 0(1). The argument in appendix (III) for the cancellation of f^= 6 - 2e from 
physical quantities works for general kinematics too. 
4 ) For n = 6, we will have surmounted the "1/e barrier" if we can produce Ie[aiajak] to 0(1), 
or alternatively the derivative d/da^ of equation (VII. 8) to 0(e). The 1^-2 t erm presents no 
problem, because we can easily compute the first derivatives of box integrals to 0(e), using 
equation (3.17) with n = 4. The I^Z^~ 2e ^ term also presents no problem, due to the manifest 
e prefactor for n = 6. Finally, the I^ =6 ~ 2e term works out as well: the 7^ term has a manifest 
e, and the rjij term requires us to know the first derivatives of I f P =6_2e to 0(e); which we can 
again compute using equation (3.17), this time with n = 6 and e — > e — 1, in terms of the 
integrals / 5 D = 6 " 2e M an d I 6 D=6_2£ through 0(1). 

There is one last step to surmounting the "1 /e barrier" for the hexagon, which is showing that 
Iq~ 6 ~ 2<e and I§ ~ 6 2e ^ drop out of "physical" quantities. Before looking at the three-parameter 
expression, let's look at the two-parameter expression (VII. 8) again and see how how Iq~ 6 ~ 2€ drops 
out of Feynman-parametrized loop integrals. Feynman parametrization of the integral l n \p^p v \ 
leads to 

1 n 

-fywz =6 - 2e m + E rf-itf-i'nM , (vii.9) 

which means that we want to show that 

E pUpu 

Because of the factor of Ag in the denominator, we should really be slightly more careful about 
how we go to "-D = 4 kinematics" , than in the one-parameter discussion above. We choose four of 
the vectors pf_i to lie in D = 4 and therefore to span D = 4; we will permit the remaining two 
vectors to have components in the [— 2e] directions, and we will only take e — ► at the end. In 
order to prove that equation (VII. 10) holds, it suffices to contract it with p^,_ 1 Pj>_ 1 , where i' and 
f each run over the set of four vectors spanning D = 4. (In the expression (VII.9) we can consider 
\x and v to belong to D = 4, not [— 2e], since we intend to contract the result with D = 4 vectors.) 
The derivation of equation (VII. 3) continues to be valid, since we are taking Pj-i to be one of the 
momenta that lie in D = 4. Thus each factor of 7, in (VII. 10) will end up with a factor of A 6 , and 
the 7^7^ term in the equation drops out in the limit e — > 0. The 77^ term has a smooth limit; using 
equation (III. 7) it is easy to show that it gives the desired result, \pv-\ 'Pj'-i- 

We now sketch how / f ? =6_2e drops out of the integral Ielp^p^p^, which after Feynman- 
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Vij 
2N f 



+ 



NrA, 



fi J 



a i a 3 = 2 5 [A}+°( € ) 



(VII.10) 



parametrization becomes 

6 1 6 

E Pf-i^-iPti / 6ka J a fe ]-2(^-2e]E^-i / 6 D=6 " 2e k] + cyclic) . (Villi) 

i,j,k=2 fc=2 

In order that the coefficient of Iq~ 6 ~ 2€ vanishes from the combination (VII. 11) we find, after 
differentiating (VII. 8), that we must have 

6 / \ , 6 



V a v X ( Wlk + Viklj + Vjkli 
/ J Pi I / /. . 1 



i,j,k=2 



2N 6 A 6 



aiOijOtk 



k=2 \^6/ 



«fc + cyclic. 



(VII.12) 

We have already dropped terms with more j^s in their numerators than Ag's in their denomina- 
tors, which will vanish in the D = 4 limit, following the same logic used earlier. If we now use 
equation (VII. 10), and again drop terms vanishing in the D = 4 limit, then we can see that equa- 
tion (VII.12) is indeed true, and so I ( f )=6_2e drops out of a momentum-space integral with three 
loop-momentum insertions. Similar considerations apply to I§~ 6 2e ^ . 

For the case of heptagon (n = 7) tensor integrals, here we will be content to obtain a well- 
defined one-parameter equation. By similar manipulations to those giving equation (VI. 9), we can 
get the one-parameter equation, 



In Wi 



,(*) 



7 



(fc) 



— 4[afc] = E-W*-i + (ji-h + 2e) 



?D=6-2e 



(VII.13) 



This equation is not adequate as it stands, since appears twice; however, by differentiating 

equation (VI. 9) with respect to a^, we get a second one-parameter equation, 

1 dL 
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n-4 + 2eda 
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(fc) 
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+ (n-5 + 2e) 
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(VII. 14) 
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(fc) 
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Solving the two equations (VII.13) and (VII. 14) for / n [aj], we get 
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(fc) 



j = l Zly n-1 
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(VII.15) 



a (k) (k) 

For n = 7 and D = 4 kinematics, we have A\> = while Nf' / 0, so we can drop the last 



term, to get: 
h[ai] 



7 (k) 



— 9/V v 

7 = 1 ZiV 6 



2V, 



(fc) 

fD=6-2£ 
(fc) J 7 



(D = 4 kinematics) . 



(VII.16) 



6 
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As with the scalar heptagon equation (VI. 10), in this equation any value of k, k = 1, . . . , 7, may 
be chosen. In using equation (VII. 16), we would like to know that I® =6 ~ 2e drops out of "physical" 
quantities. This amounts to showing that 

Erf-if^fe = °. ( VIL17 ) 

i=2 ZiV 6 

for D = 4 external kinematics. To show that equation (VII. 17) holds, we contract it with four 
independent vectors spanning D = 4 Minkowski space, namely V^-\ for j £ {2, . . . , n}, j / fc. We 
thus have to show 

7 (fc) 

ES) Pi-i-Pi-i = 0. (VIL18) 

But this is the same sum encountered in showing that D ~ 6 2e drops out of "physical" linear 
combinations of i^\ai\, which we have already shown above. 

Finally, the linear combinations of i^[ai] that appear explicitly in equation (VII. 16), namely 
Yli=2Pi-i ai a *' are a ^ S0 ^ ne same as those occurring in "physical" one parameter hexagon integrals 
(using Iq^ [a,j] = 0). So 1^ D ~ 6 2e drops out there too. Therefore "physical" combinations of Ij[ai\ 
in equation (VII. 16) are given in terms of well-defined, D = 4 quantities, as desired. 

To get heptagon integrals with two or more Feynman parameters inserted, one can differentiate 
equation (VII. 15) with respect to the CKj, and then take the limit of D = 4 kinematics; it remains 
to show that the unwanted six-dimensional integrals drop out for "physical" quantities. 
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Figure Captions 

Fig. 1: A schematic depiction of equation (3.24), with the coefficients suppressed. 
Fig. 2: A diagram containing a triangle loop with one massive (or off-shell) leg. 
Fig. 3: A diagram containing a box loop with one massive leg. 
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Fig.l 




Fig. 3 



